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ABSTRACT 

The rotation-mass-age relationship offers a promising avenue for measuring the ages of field stars, 
assuming the attendant uncertainties to this technique can be well characterized. We model stel- 
lar angular momentum evolution starting with a rotation distribution from open cluster M37. Our 
predicted rotation-mass-age relationship shows significant zero-point offsets compared to an alterna- 
tive angular momentum loss law and published gyrochronology relations. Systematic errors at the 
30 percent level are permitted by current data, highlighting the need for empirical guidance. We 
identify two fundamental sources of uncertainty that limit the precision of rotation-based ages and 
quantify their impact. Stars are born with a range of rotation rates, which leads to an age range at 
fixed rotation period. We find that the inherent ambiguity from the initial conditions is important 
for all young stars, but becomes large for old stars below 0.6 M©. Latitudinal surface differential 
rotation also introduces a minimum uncertainty into rotation period measurements and, by exten- 
sion, rotation-based ages. Both models and the data from binary star systems 61 Cyg and a Cen 
demonstrate that latitudinal differential rotation is the limiting factor for rotation-based age precision 
among old held stars, inducing uncertainties at the ~ 2 Gyr level. We also examine the relationship 
between variability amplitude, rotation period, and age. Existing ground-based surveys can detect 
field populations with ages as old as 1 - 2 Gyr, while space missions can detect stars as old as the 
Galactic disk. In comparison with other techniques for measuring the ages of lower main sequence 
stars, including geometric parallax and asteroseismology, rotation-based ages have the potential to be 
the most precise chronometer for 0.6 - 1.0 M© stars. 

Subject headings: Stars: evolution - Stars: rotation - Stars: late-type 
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1. INTRODUCTION 

Ages of stars are both fundamental and difficult to 
measure. Stars of known age, such as members of star 
clusters, have therefore been studied intensively. For the 
far more numerous field stars, by contrast, our tools for 
inferring ages are blunt and inexact. We can infer ages 
for some binary stars, and also for turnoff stars, from the 
combination of parallax and metallicity measurements. 
However, for typical M-G dwarfs, slow luminosity evolu- 
tion on the main sequence and ambiguities from metal- 
licity and helium errors result in uncertain age estimates. 
Rota tion provides an attrac t ive alternativ e chronometer 
(e.g. iLachaume et al.l 119991: iBarnesI 120071) . Stars with 
minimal nuclear evolution can experience strong rota- 
tional evolution, in the sense that low mass stars are 
observed to rotate more slowly as they age. The idea of 
rotation as a clock for low mass stars is an old one; how- 
ever, the prospect of rotation as a precise clock is new. 
We identify two effects that limit the intrinsic precision of 
rotation-age relationships: the large initial range of stel- 
lar rotation rates and surface differential rotation. We 
dedicate this paper to a comparison of rotation-mass- 
age predictions and a critical evaluation of the limiting 
precision of the rotation-mass- age relationship. 

1.1. A Standard Picture of Angular Momentum Loss 
in Low Mass Stars 
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Measurements of Ca H and K emission for main se- 
quence stars led [Wilson (1963) to insightfully notice an 
inverse relationship between chromos pheric activ i ty and 
stellar age. In an influential paper, ISkumanichl ()1972f ) 
argued that stellar rotation, activity, and lithium abun- 
dances for solar-like stars all decrease as the square root 
of age. Further progress was hindered by basic observa- 
tional and theoretical limitations. 

The theory of angular momentum evolution is complex. 
To predict rotation-age relationships, one must under- 
stand the origin of stellar rotation, angular momentum 
transport in stellar interiors, and angular momentum loss 
from magnetized solar-like winds. From an observational 
point of view, it is difficult to measure rotation periods in 
old field stars from either photometric or spectroscopic 
methods. This difficulty has been amplified by the even 
smaller subset of stars with independent age estimates. 
A brief summary of the lessons learned is useful for plac- 
ing ages inferred from rotation into the modern context. 

Open clusters and star forming regions have been cru- 
cial for developing the standard picture of angular mo- 
mentum evolution in low mass stars. Large rotation 
datasets, and models capable of a ddressing them, be- 
came available in the la te 1980s (Stauffer et al. 1989; 
iPinsonneault et aLlll989f ). T Tauri stars represent the 
initial hydrostatic evolutionary stage for stars, and rota- 
tion measurements in star forming regions therefore set 
the initial conditions for angular momentum evolution. 
Open clusters provide snapshots in time of the evolution 
of stellar rotation rates. The simplest model of angular 
momentum evolution would begin with the distribution 
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of rotation rates in protostars and evolve them forward 
assuming solid body rotation and a scaled solar angular 
momentum los s rate dJ/dt ~ w 3 (|Kawalerl 1988), which 
reproduces the iSkumanichl law. Nature is not obliged to 
conform to such an idealized formula, however, and the 
obstacles to a more complete theoretical description are 
severe. Empirical guidance is thus essential. 

A standard working model of angular momentum evo- 
lution has emerged, with some important refinements on 
the simplest possible model. Protostars of the same mass 
are b orn with a range of rotation rate s (e.g.lBouvier et all 
1986b for more recent w ork see [Rcbull c t al.l 12004 
Henderson fc Stassunl 12011ft. Early models starti ng in 
the protostellar regime ()Pinsonneault et all 119901 ) pre- 
dicted that stars would converge to a narrow mass- 
rotation relationship at fixed age, but did not correctly 
predict the existence of either the most rapid or the slow- 
est rotating stars in young open clusters. Three impor- 
tant refinements were introduced to explain the observed 
angular momentum evolution pattern. 

The survival of rapid rotation in young stars prohibits 
a simple scaling of the solar torque by the cube of the ro- 
tation rate; saturation in the loss rate is required. A sim- 
ilar behavior is observed in proxies of ma gnetic activity 
(jNoves et al .l 119841: iPizzolato et all 12003ft an d expected 
in wind models (JMacGregor fc Brennerlll991l ). Both the 
activity and spin do wn data require a mass dependent 
saturation threshold (jSills et al. 2000), in the sense that 
lower mass stars saturate at lower rotation rates. 

The opposite problem occurs for the slowest rotation 
rates seen in young open clusters; scaled solar winds are 
too inefficient to reach such low rotation r ates so quickly . 
This led to the h ypothesis of disk-locking (|Koenigllll991l : 
iShu et all 119941 ) . which is the hypothesized rotational 
coupling between the stars convective envelope and the 
inner accretion disk. Observations require the existence 
of a net torque in excess of that expected from solar- 
like winds, and that this effect is more important for 
some stars than others. Characteristic lifetimes of order 
6 Myr and les s are observational ly plausible f or gaseous 
disk lifetimes (Hai scti et alll200ll but see also iMamajeld 
2009) , and they can re concile the protoste llar and Z AMS 
rotation distributions ((Tinker et al.ll2002ft . 

The final observational effect is the detection of core- 
envelope decoupling. Stars spin up as they contract onto 
the main sequence, where their structure stabilizes and 
the timescale for evolution increases. A magnetized wind 
will initially spin down the convective envelope; the ra- 
diative core will then respond over the internal angular 
momentum transport timescale. There will be a rapid 
initial drop in surface rotation if the transport timescale 
is long enough, followed by an epoch where the core 
spins down and the envelope does not (e.g. lEndal fc Sofia! 
1979). Once the entire star is coupled it will spin down as 
a solid bod y, even if there is a differen tially rotating inter- 
nal profile (|Pinsonneault et"aTlll9 89). The observational 
signature is that the slow rotator population in young 
open clusters spins down quickly at first, and then the 
spin down sta lls. Numerous studies have confirm ed this 
phenomenon (jK cppcn set al.lll99 5: Krishnamurt hi et al. 
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spin down as a solid body, suggesting that the transport 
timescale is a function o f both rotation rate and mass 
(|Denissenkov et all 12010ft . The underlying mechanisms 
for all three of these effects are active research topics, 
but they are all physically well-motivated and supported 
by strong empirical evidence. A practical consequence 
is that there are documented problems with adopting a 
Skumanich-style spin down model in young stars. 

1.2. Rotation Data and Gyrochronology 

The complex physical picture for young stars is daunt- 
ing. Fortunately, there has been an explosion in the 
quality and quantity of stellar rotation data. This 
has largely been driven by planet transit searches, but 
dedicated stellar rotation programs, such as the Mon- 
itor P ro ject, ha ve al so been valuable (llrwin e t al. 2006, 
l2m7llb1 . l2OT8lbl and lMeTbom et aLII2009LI2011allbft . Stel- 
lar rotation periods are now available for bulk stel- 
lar populations, especially for cool and rapidly rotat- 
ing stars. Space missions have extended our precision 
in both amplitude and period, raising the prospect of 
accurate rotation periods, even for old Sun-like stars 
(|Basri et al.l 12010ft . The qualitative rotation-mass-age 
relationships from earlier decades are now being re- 
evaluated as practical tools for stellar population stud- 
ies, and open cluster stars have once again become a 
central calibrating tool. Modern rotation-mass-age, or 
gyrochronology, relations (IPace fc Pasauin i 2004; Barnes 
poUlMamaiek fc Hillenbrandll2008l;IMeibom et al.ll2009L 
l2011al) are based on the following major insights gleaned 
from open cluster studies. 

Young open clusters reveal that coeval stars of the 
same mass are born with a range of rot ation rates (e.g. 
NGC 65 30. Henderson fc Stassunl I20TTI: Orion Nebula 
Clust er, lAttridge fc Herbstl 11993 IC 348, IHerbst et all 
12000ft . However, the majority of young stars arrive on the 
main sequence as slo w rotators (Stauffcr et al.lll989l and 
lHartman et alll2010l are nice examples for the Pleiades 
cluster.) One can therefore define a median relationship, 
even in young systems, which includes many of the stars. 
The width of the rotation distribution also narrows with 
time; this process is fastest for solar analogs and takes 
progressively longer for lower mass stars. This conver- 
gence is seen in solar analogs by ages of order 500 - 600 
Myr, exemplified by the Hy ades (Radick ct al. 1983) and 
M37 partman et all 12009ft . Finally, the median rota- 
tion period at fixed color increases with increased age. 
Published gyrochronology prescriptions take advantage 
of these basic features to define working relationships 
between rotation, color, and age using a semi-empirical 
approach. The median rotation of the bulk population 
is defined for young open clusters, tested in clusters with 
ages under 1 Gyr, and extrapolated to older st a rs us- 
ing the solar calibration as a constraint. iBarnesI (|2007ft 
suggested that the influence of color (mass) and age are 
independent effects and thus mathematically separable 
into two functions as shown below: 



2010ft . with inferred timescales in the range of tens to 
hundreds of Myr, depending on mass. The upper en- 
velope of the rotation distribution appears to always 



P(B-V,t) = f(B-V)g(t) 
f(B-V) = a[(B-V) -c}' 



(i) 

(2) 
(3) 



where the age t is in Myr. Spin down is accomplished 
with a power law g(t) and an index of n = 0.5 corre- 
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TABLE 1 

Parameters in Gyrochronology Models 



Source 



Barnes (2007^ 
Mamajck & Hillenbrand 
Meibom et al. (2009 ~ 



0.7725 ± 0.011 
0.407 ± 0.021 
0.770 ± 0.005 



0.601 ± 0.024 
0.325 ± 0.024 
0.472 ± 0.027 



0.40 

0.495 

0.553 



0.010 
0.052 



0.5189 ± 0.0070 
0.566 ± 0.008 
0.52 



sponds to the iSkumanichl (|1972l ) law. The polynomial 
/(£) models the color dependence of the slow rotator se- 
quence at fixed age. Since f(B — V) lacks time depen- 
dence, any relationship of this form produces self-similar 
mass-rotation isochrones. Although many variations ex- 
ist, for the purposes of this paper, we compare three 
published gyrochronology relationships, with parame- 
ters listed in Table [TJ These gyrochronology relation- 
ships are typically calibrated to clust er data i n the color 
range 0.5 < (B - V) n < 1.4. Both iBarnesI ()2007f ) and 
iMamaieTT fe Hillenbrand (2008) base the color-rotation 
relationship on many well studied clusters (e.g. a Per, 
the Pleiades, Hyades, etc.) and determine the index n 
as the factor t hat reproduces the Su n's rotation rate at 
the solar color. iMeibom et al.l J2009) fit the color depen- 
dence in a single cluster ( M35) , but retained the power 
law index n from Barnes (2007). 

Gyrochronology is a promising concept, but there is 
real tension between our understanding of angular mo- 
mentum evolution, cluster data, and a simple mapping 
of rotation onto a mass-dependent age. This is especially 
true for rotation-mass- age relationships in the context of 
young clusters, which we examine in a companion pa- 
per (Epstein, Pinsonneault, & Terndrup 2012; Paper II.) 
Fortunately, there is strong evidence that the physical 
picture is simpler for the older stars that will dominate 
a field sample. Protostar-accretion disk interactions are 
important on timescales of millions of years, while ef- 
fective core-enve lope coupling is establis hed by an age 
of 500 Myr (e.g. IDenissenkov et al.|[20i0l) . The bulk of 
the well-studied open cluster systems are clustered in age 
around the Pleiades (125 Myr) and Hyades (625 Myr), 
with a few systems in between. 

In this paper, we therefore focus on the expected be- 
havior of old field stars, with an emphasis on stars older 
than 500 Myr. The distribution of initial conditions can 
be empirically specified by open cluster stars, and the as- 
sumption of solid body spin down at this age and beyond 
is well-justified. Furthermore, there are good constraints 
on angular momentum loss, permitting a well-posed for- 
ward modeling of angular momentum evolution. 

Much recent work has focused on presenting gy- 
rochronology recipes as a function of color. We compare 
these relationships with angular momentum loss mod- 
els and the extremely limited current data for old stars. 
New data, especially from the Kepler mission, will shed 
valuable light on the actual functional form of rotational 
evolution in old stars (a point that we will return to in the 
conclusions.) Here, we emphasize an under- appreciated 
aspect of the problem: how precise are rotational ages? 
For the purposes of this question, we assume that the 
mean rotation-mass-age relationship may be sufficiently 
well calibrated in the future such that its uncertainty 
is not the limiting factor in the error budget. Instead, 



we investigate the irreducible sources of uncertainty in 
rotation-based age measurements. In particular, both 
the initial scatter in stellar rotation and error in rota- 
tion period measurements place fundamental limits on 
rotation-based age measurements for old field stars. 

This paper is organized as follows. Section [5] contains 
a detailed description of the initial conditions, stellar in- 
terior model, braking law, and assumptions necessary to 
model stellar angular momentum evolution. We present 
our results in Sj3j We compare the median rotation- 
mass-age relationship predicted by our preferred angu- 
lar momentum loss model with an alternate braking law 
( t|3.1.1[) and published gyrochronology relations ft j3.1.2p . 
Those theoretical predictions are tested against observa- 
tional constraints from two binary star systems 61 Cyg 
and a Cen in ^3.1.31 Sources of uncertainty in obtaining 
rotation-based ages are quantified in ^3.21 These include 
both fundamental limitations (e.g. arising from the initial 
scatter in rotation rates, N3.2.H and from surface differ- 
ential rotation, i )3.2.2p and model-specific errors ( i)3.2.3l 

MM- 

Kepler and CoRoT are measuring asteroseismic ages 
for tens of thousands of old solar- like stars (e.g. 
IMosser et allllHOt [GlhTland et al.|[20Tol ). However, such 
space missions cover only a limited number of fields and 
targets. We therefore also examine how the prospects for 
rotation-mass-age relationships depend on the sensitivity 
of the survey. In i j3.3[ we demonstrate that ground-based 
studies can find young to intermediate aged populations 
(based on variability) for plausible precision ranges, while 
it is likely that the detection of older populations will re- 
quire a degree of precision currently achievable only from 
space. Finally, in 21 we compare the limiting precision of 
rotation-based ages with other age diagnostics: trigono- 
metric parallaxes and asteroseismology. 

2. METHOD 

Because of the paucity of data at old ages, we rely on 
the existing theoretical framework to predict the behav- 
ior of the rotation-mass-age relationship at late times. 
As described in fll.H observations of rotation in open 
clusters (< 1 Gyr old) have guided theory to a standard 
model of angular momentum evolution. We treat the 
time evolution of rotation as an initial value problem. 
The initial mass-rotation distribution for our models is 
empirically determined by meas urements of the po pulous 
open cluster M37 (550 Myr; lHartman et al.ll2009l ) Stel- 
lar angular momentum evolution is modeled with solid 
body angular momentum transport and a braking law 
specifying the torque provided by a stellar wind acting 
on a star. We compare two prescriptions for the strength 
of the wind: the commonly-use d modified Kawaler law 
(jKawalerll 19881: ISills et al.ll2000f) and the re cently formu- 
lated loss law by [Reinersl&l^oha^a (HHQ) ■ These rota- 
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tional evolution models piggyback on th e underlying Yale 
Rota t ing Evolutionar y Code (YREC; iDemarque et al.l 
120081 iSills et aT1l2000D whose interior models determine 
the evolving stellar structure. 

2.1. Initial Conditions 

One could, in principle, use the rotation distribution 
of T Tauri stars and model protostellar-disk interactions 
to obtain the initial conditions. Such an exercise is valu- 
able when testing the origin of the distributions seen in 
young systems or when studying core-envelope decou- 
pling. However, our focus is to investigate the impact of 
the observed spread of rotation rates on age estimates. 
We bypass model-induced uncertainties in the initial ro- 
tation distribution by adopting an empirical rotation dis- 
tribution drawn from measurements of stellar rotation in 
a cluster. 

Selecting the age to which models are initialized re- 
quires careful consideration. Stars do eventually forget 
their initial conditions as they spin down and converge to 
a common rotation-age sequence. Among clusters older 
than ~ 500 Myr, the mass-rotation sequences for differ- 
ent clusters are consistent with being drawn from the 
same parent population (see Paper II). This commonal- 
ity suggests that results obtained for a cluster older than 
500 Myr may be extended to both other clusters in the 
same age regime and field stars. Using empirical distri- 
butions to set the initial conditions, however, does mean 
that the results are sensitive to the same biases and sys- 
tematics present in the input cluster data. To mitigate 
this problem, we use a rich cluster with well understood 
biases, measures of completeness, and little field con- 
tamination, namely the ope n cluster M37 (NGC 2099). 
lHartman etafl (J2008, 2009) performed a deep, 20 night 
survey with the MMT of M37, reaching down to r ~ 23 
mag. They determined an age of 550±30 Myr (with over- 
shooting), metallicity [M/H]=+0.045 ± 0.044, and dis- 
tance modulus (m—M)v — 11.572±0.13. reaching down 
to r ~ 23 mag. Searching for pe riodic va riabilit y using 
a mu ltiharmonic AoV algorithm (Schwarzcnbcrg-Czerny 
1996) yielded a clean sample of 372 rotation period mea- 
surements. The sample is close to 100% comp lete to 
r = 20 mag (see Figure 9 in lHartman et al.ll2009|) . in the 
sense of detecting stars with modulation above 0.01 mag, 
although the color-period diagram may be biased toward 
brighter stars. This means that an increasing fraction of 
cluster members have undetectable variations as mass 
declines. Hartman et al. estimate ~ 20% contamination 
o f the cluster sam ple by background field stars. 

lHartman et al.l obtained B, V, and Ic band photome- 
try for this cluster. For a given set of isochrones, dif- 
ferent colors resulted in consistent mass estimates in 
regimes where the isochrones were not double-valued. 
For this paper, we select (V-Ic)o to estimate mass be- 
cause it has a large dynamic range and its relationship 
to mass remains monotonic for cool stars (see Figure [T]). 
The (V-Ic) photom etry has been correc ted for redden- 
ing E(V-Ic)=0.355 (|Hartman et al. 2008) and compared 
to the corresponding 550 Myr, [Fe/H] = +0.045 YREC 
isochrone. 

As demonstrated by Figure [TJ the choice of isochrone 
induces larger shifts to the mass scale than the choice of 
color. The color-mass relationship for different isochrone 
systems agrees for bluer stars (masses > 0.75 M Q ), 




Fig. 1. — Color-magnitude diagrams (top) of M37 with YREC 
(solid) and Padova (dashed) [Fe/H] =+0.045 isochrones at 550 Myr 
overlaying the cluster main sequence (gray points). The corre- 
sponding color-mass relationships appear in the lower panels. 

but breaks down in the red. The YREC and Padova 
isochrones were selected because their differences pro- 
duce the largest shift in the mass scale between any two 
widely-used isochrone systems. We note that the YREC 
isochrones better reproduce the morphology of the clus- 
ter main sequence compared to the Padova isochrones, 
so our error estimates will be conservative. 

2.2. Interiors 

We anchor our angular momentum study in theoretical 
models of stellar structure calculated using the YREC 
stellar interior models. We compute solar metallicity 
models (X=0.7190, Z=0.0165) with no diffusion. These 
provide a prior on a star's He abundance, radius, moment 
of inertia, convective overturn timescale (tcz), and lu- 
minosity as a function of stellar mass and age. 

The standard model of angular momentum loss de- 
scribed in Ml . II is only valid for stars in a particular mass 
range. Because stars above spectral type F8 show no 
obvio us sign of a rotation-age relationship (IWolff et all 
119861 ). this sets an upper mass limit on where rotation- 
based ages are well motivated. The old open clusters 
have sparse data at the low mass end, and there is 
evidence for a departure from a R ossby scaled sat ura- 
tion threshold in that mass regime (|Sills et al.l 120001 : see 
i )2.4.1[) . Furthermore, quantities like tcz become ill- 
defined for fully convective stars. We therefore restrict 
our calculations to the mass range of 0.4 M Q < M < 
1.2 M , sampled in equally spaced 0.05 M Q intervals. 
Stars in M37 (used to define the initial rotation distri- 
bution, ij2.ip are excluded if their mass lies outside of 
this restricted range where our angular momentum loss 
model is valid. 

The empirically determined stellar masses and rota- 
tion rates from M37 set the initial conditions ( §2.1|) . A 
prescription for angular momentum loss ( H2.4p is used 
to simulate M37's rotational future. We calculate a 
star's angular momentum at each point in time on a 
logarithmically-spaced grid of ages. At fixed age, we per- 
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form a polynomial interpolation in logarithmic rotation 
period on a grid of models to the mass of the star. 

We set the birth line by requiring that the deuterium 
abundance declines to a level of 0.1% below the initial 
value, resulting in a slight offset from the iSills et all 
(2000) models. We stop our models if their age exceeds 
either a stars' lifetime on the main sequence or the age 
of the Galactic disk, assumed to be 10 Gyr. The main 
sequence turnoff is defined as the point when the core 
hydrogen abundance drops below Xc = 10~ 3 . 

In accordance with standard mixing length theory, we 
define tqz locally as tcz — H p /v conv evaluated one pres- 
sure scale height above the base of the surface convection 
zone. In the YREC interior calculations, there is numer- 
ical jitter in tqz at the < 1% level over a stars's main 
sequence lifetime. We therefore average tqz over the 
five closest points in time to smooth over these numeri- 
cal perturbations in the depth of the convection zone. 

2.3. Internal Angular Momentum Transport 

Core and envelope coupling timescales are short in 
comparison to th e lifetime of G, K, and M dwarfs dis- 



cusse d here (e.g. iDenissenkov et all 1201 0; Spad a~et~aTl 
1201 lh . Because we are studying the rotation of field 
stars, assuming that the coupling timescale for angular 
momentum transport between the core and envelope is 
short (i.e. solid body rotation is enforced at all times) 
is a reasonable approximation for the majority of the 
age range of interest. This well motivated simplifica- 
tion has the added benefit of restricting the number of 
free parameters in the angular momentum loss law. The 
rotation rates in young systems are kno wn to be influ- 
enced by core-envelope decoup l ing (e.g . Kcppcn s et al 
1995; Kri shnamurthi et al.lfl 998; Allain 1 9981: llrwin et al 
2007b; iCollier Cameron et al.l l2009t IDenissenkov et al. 
2010f ) and one purpose of our work is to infer the impact 



on global angular momentum evolution studies. We dis- 
cuss tests of the coupling timescale by comparing results 
from clusters with ages close to and much older than the 
t c in Paper II. M37, the cluster used here to set initial 
conditions, falls into the latter category. 

2.4. Torque 
2.4.1. Kawaler Wind Law 

To simulate stellar angular momentum loss, we adopt 
as our preferred model a modified K a waler wind law 
(jChabover et al.lll995l ISills et alll2000h IKawalerl (fl988h 
was able to recover the iSkumanichl (|1972f ) law (cj ex 
i" 1 / 2 ) for his favored choice of the magnetic field ge- 
ometry. This yields a wind law of the form dJ/dt = 
—Kuj d 'R 1 /' 2 M~ 1 /' 2 1 where the normalizing constant K 
parametrizes our ignorance about winds and magnetic 
field generation in stars. 

As initially conceived, this law does not account for 
magnetic saturation- a necessary element in the stan- 
dard theory of an gular momentum loss ( 31-1[) . We follow 
ISills et all (J2000I ) in modifying the original Kawaler loss 
law to include a prescription for how the torque changes 
when a star's magnetic field strengt h saturates at a crit- 
ical ro tation rate u> cr it- Because the lReiners fc Mohantvl 
(2012) wind law (discussed in the next section N2.4.2p 
treats magnetic saturation differently, we supply a brief 
discussion of the empirical motivation and prior theo- 



retical work that produced the functional form of the 
saturation threshold adopted here. 

Observational studies of rotation and magnetic activity 
indicators support a Rossby-scaled magnetic saturation 
threshold. Both chromospheric Ca II H and K and coro- 
nal X-ray emission display a break in their relationship 
with Ross by number, where the Rossby number is R g = 
Prot/rcz (jNoves et al.lll984l: iPatten fc Simonlll996D . In 
fact, activity correlates m ore tightly with Ro ssby num- 
ber than rotation period (jRadick et al.l 11987) . Because 
measures of stellar activity serve as a proxy for magnetic 
field strength, which in turn powers the wind, the rate of 
angular m omentum loss is also though t to plateau. Con- 
sequently, Krishna murthi et"al"1 (|1997| ) proposed the ad- 
dition of a Rossby scaled saturation threshold for stars 
more massive than 0.5 M Q , resulting in the following 
modified Kawaler wind law: 




-•^kU< 






where 



' J crit — Ldcrit,Q~ 



T CZ,Q 
TCZ 



(5) 



From the solar mass YREC model, tcz,q = 1-01 x 10 6 
s. Because lower mass stars have deeper convection 
zones, Tcz contains the mass dependence in the satu- 
ration threshold. 

X-ray measurements confine t he saturation threshol d 
QJerit to the ran ge 5 to 20 w Q (|Patten fc Simonl 1 1996ft . 
ISills et all (|2000D tested the Rossby scaling by comparing 
theoretical models to rotation distributions of IC 2602, 
IC 2391, a Per, the Pleiades, and H yades. They found 
that the iKrishnamurthi et all (|T997) scaling (uJ cr it,Q = 
10 ujq = 2.86 x 10 -5 s _1 for solid body models) works 
for stars with masses between 0.6 M Q < M < 1.1 M Q . 
Howev er, the Ros sby scaling breaks down for lower mass 
stars; ISills et all argue that the saturation threshold 
drops more steeply than a Rossby scaling for stars with 
0.1 M < M < 0.5 M . They then used the clus- 
ters to empirically measure the mass dependence of u! cr it 
for stars in this mass range. For solid body models, 
they find that us cr it t Q = 7, 6.2, 5.1, 3.6, and 1.7 ujq at 
0.5, 0.4, .3, 0.2, and 0.1 M Q . We interpolate on the 
ISills et al| calibration for stars M < 0.5 M Q and adopt 
the IKrishnamurthi et al.l calibration for M > 0.5 M Q . 

Given this definition of ui cr it , Eq. U is calibrated such 
that a 1 M Q star achieves the solar rotation rate ojq = 
2.86 x 10~ 6 rad s _1 (the Sun's equatorial rotation period 
is P = 25.4 days) at the solar age of 4.5 Gyr starting 
with an initial period of 10 days and no disk locking. 
This results in K Q = 3.15 x 10 47 g cm 2 s. The choice 
of initial rotation period is plausible; it is similar to the 
av erage T Tauri rotatio n period of 9.45 days measured 
bv lChoi fc Herbstl (J1996D in the ONC. However, the pre- 
dicted rotation period at the age of the Sun, and there- 
fore the value of the normalization constant, is highly 
insensitive to the choice of initial rotation period. This 
is a characteristic feature of angular momentum laws of 
with a strong dependence on rotation rate, namely, that 
once stars enter the unsaturated regime (w < w cr it), they 
rapidly forget their initial conditions. For example, if 
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Fig. 2. — Rotation period evolution according to the modified 
Kawaler (Eq. [4] solid) and Rcincrs & Mohanty (Eq. [6] dashed) 
braking laws for a 0.5 Mq (top) and 1 Mq (bottom) star. All 
models assume an initial rotation period of 8.5 days. Evolutionary 
trac ks are shown for scen arios with and 5 Myr of disk-locking. 
The IReiners & Mohanty| law has been re-calibrated such that it 
reproduces the Sun with the YREC interior models. The diagonal 
dotted line gives the solar calibrated Skumanich P ~ t 0,5 relation, 
calibrated to the Kawaler law at 550 Myr with 5 Myr disk-locking 
(top) and to the Sun (bottom). The horizontal dotted line shows 
the mass- and time- independent saturation threshold adopted by 
Rcincrs & Mohanty (2012). Grey solid lines show the case of no 
braking for comparison. 



we instead calibrated the solar rotation rate with an ini- 
tial period of 10 days and a 5 Myr disk locking time, this 
would result in a slightly different normalization constant 
K Q = 2.97 x 10 47 g cm 2 s. 



2.4.2. \Reiners & Mohantil Wind Law 

iReiners &: Mohantvl (|2012f ) present a new braking law 
that pivots on two major deviations fro m the modified 
Kawaler law. First. IReiners fc Mohantvl argue that rota- 
tion is related to the magnetic field strength (Bo oc us a ) 
rather than the surface magnetic flux (BqR 2 k w°) as 
assumed in [Kawaler (19881). This factor of R 2 difference 
has powerful ramifications when propagated to the an gu- 
lar momentum loss law. Second, IReiners fc Mohantvl as- 
sume that the saturation threshold uj cr n does not depend 
on mass. This differs from the Rossby-scaled saturation 
threshold used in §2.4.1) in this new law, the saturation 
threshold is insensitive to structural differences in stellar 
interio rs. The reader should refer to IReiners &: Mohantvl 
(2012) to see the specifics of their derivation, but the 
resulting angular momentum loss law (their Eq. 5) is re- 
produced below: 



dJ 
~dt 




for w > u; crit 
for uj < w crit 



(G) 



where C is a constant relating to the magnetic field 
strength and rate of mass loss, R is stellar radius and M 
is stellar mass in cg s units. It is important to note that 
IReiners fc Mohanty use the saturation threshold in the 



opposite sense of the modified Kawaler law. In the satu- 
rated regime (to > aj cr it), the rate of spin down is set by 
a maximum magnetic field strength (B cr ; t ) for all stars. 
The saturation threshold 0J cr it only affects the torque in 
the unsaturated regime (uj < uj cr u). This means that 
changing the saturation thr eshold affects the la t e-time 
behavior of stars under the IReiners fc Mohantvl (|2012f ) 
prescription. This is not true for a modified Kawaler 
law, where the effect of u; cr it is only felt in the saturated 
regime and erased as stars forget the ir initial co nditions. 

To calibrate this wind law, [Reiners &: Mohantvl match 
the solar rotation rate at the solar age, starting with 
an initial rotation rate of 8.5 days. The presence of a 
disk is assumed to hold the rotation rate constant for 5 
Myr. This yields their best fit combination of parameters 
[ Wcrit , C] = [8.56xl0- 6 s" 1 , 2.66xl0 3 (gVcm^s 3 ) 1 / 3 ]. 
We recalibrate the constant C to reproduce the Sun with 
our models and a 5 Myr disk-locking lifetime, finding a 
different value of C = 3.16xl0 2 (g^cni^s 3 ) 1 / 3 . 

Figure [2] compares angular momentum evolution un- 
der both a modified Kawaler and a Rei ners fc Mohantvl 
wind law with and without a disk. The largest differ- 
ences ap_p_ear_among_Jow mass stars. As discussed in 
§2~4~TI iSills et all (|2000D find that t he Rossby scaling is 
inadequate for stars < 0.6 M©. The IReiners fc Mo hanty 
law naturally preserves rapid rotation of low mass stars. 
For solar mass stars, the Kawaler wind law asymptotes 
to the Skumanich relationship, becoming nearly indistin- 
gui shable by 500 Myr in the disk-locked case. In contrast, 
the IReiners fc Mohantvl wind law predicts a slower spin 
down than Kawaler does for the first few Myr, followed 
by an interval of rapid angular momentum loss between 
~ 200 — 700 Myr, before leveling off and approaching a 
rate slightly shallower than Skumanich later in life. 

3. RESULTS 

Because the timescale for rotational evolution is faster 
than that of nuclear burning for low mass stars, rotation 
represents a promising tool for measuring stellar ages. 
Here we explore the potential of using rotation as a clock 
for field stars. Having initialized our clock in §2J we now 
test how it runs compared to other timepieces. Section 
l3.1l compares our models' median rotation-mass-age rela- 
tionship with published gyrochonology relations. These 
theoretical predictions are then evaluated by their agree- 
ment with the independently-measured age of two binary 
star systems, and also with slowly rotating field stars. 
Such empirical guidance is necessary for discriminating 
between different theoretical prescriptions. 

Until now, clockmakers have focused on how well 
rotation-based clocks keep time, but largely ignored the 
precision of those measurements. In §3.21 we investi- 
gate many sources of error that can affect rotation-based 
ages. The dominant source of uncertainty for low mass 
stars proves to be the intrinsic scatter in the initial ro- 
tation distribution. Because these stars require many 
Gyr to converge to a unique rotation-mass-age relation- 
ship, this induces large uncertainties in rotation-based 
ages. Among old stars, the range of surface rotation rates 
observable due to latitudinal differential rotation limits 
the precision. Systematic uncertainties due to modeling 
choices are also discussed. 

Stars, like mechanical clocks, wind down over time. 
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More slowly spinning stars tend to be less magneti- 
cally active and t heir amplitude of variability declines 
(|Skumanichl I1972D . Past a certain age, the variability 
will become too small to be detectable from the ground. 
In A3.31 we determine that this technique of measuring 
rotation-based ages is feasible for stars as old as 1-2 Gyr 
using groundbased surveys, but space-based instruments 
are required to see photometric modulation in older pop- 
ulations. Comparison of our predictions with rotation 
periods detected by current groundbased surveys sug- 
gests that our predictions may be too conservative and 
that such studies may probe even older populations. 

3.1. Median Rotation-Mass- Age Relationships 

Because currently-studied open clusters only provide 
guidance for stellar angular momentum evolution < 1 
Gyr, we utilize the angular momentum loss models de- 
scribed in $2] to forecast rotational behavior over the life- 
time of the Galactic disk. Section 13.1.11 characterizes 
how both the median and shape of the rotation distri- 
bution evolve over time for two angular momentum loss 
models. The rate at which stars converge to a unique 
rotation-mass-age sequence sets the timescale for when 
rotation becomes a useful clock. We compare the median 
rotation-mass-age trends for these models with published 
gyrochronology relationships in M3.1.2I Because some of 
the gyrochronology relationships fail to reproduce the 
median trend seen in the M37 cluster data, we question 
the calibration of these relationships. 

Although observational constraints on the rotation- 
mass-age relationship are rare at late times, two binary 
star systems, 61 Cyg and a Cen, possess independent age 
measurements. We therefore test the rotation-based age 
predictions from our models and published gyrochronol- 
ogy relationships against these binaries in §3.1.31 In the 
case of these old binary stars, latitudinal surface differ- 
ential rotation induces ~ 2 Gyr uncertainties in rotation- 
based age estimates. We also consider the implications 
for these angular momentum evolution models of long 
rotation periods observed for old field stars. 

3.1.1. Time Evolution of M37's Rotation Distribution 

Young clusters display a large spread in stellar rota- 
tion rates across all masses. This wide range converges 
as stars age, but the timescale for this process can be 
long, especially in lower mass dwarfs. We begin by mod- 
eling the time evolution of M37's rotation distribution 
and investigating how long it takes for rotation rates to 
converge. The timescale for convergence is important 
because the median rotation-mass-age relationship only 
becomes a good chronometer after stars have forgotten 
their initial conditions. 

The description in §2] provides the theoretical frame- 
work necessary to predict the rotation of a star at any 
future time, given its mass and initial period. We use 
solid body models to evolve M37's rotation distribution 
from 550 Myr forward in time. We explore two differ- 
ent prescriptions for the torque, starting with a modified 
Kawaler wind law ( §2.4.ip . Figure [3] shows M37's initial 
rotation distribution at 550 Myr and also how that dis- 
tribution would evolve under the influence of a modified 
Kawaler wind law. For visual clarity, we show the pro- 
jected distribution at representative ages 1, 2, 5, and 10 
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Fig. 3. — M37's color-period distribution (Hartman ct al. 2009) 
evolved forward in time from 0.55 Gyr (red stars) to 1 Gyr (gold 
squares), 2 Gyr (green triangles), 5 Gyr (blue circles), and 10 Gyr 
(purple crosses). The corresponding ages (in Gyr) are indicated 
to the left of the mass-rotation sequence. These synthetic rotation 
distributions are created by evolving the M37 rotation distribu- 
tion forward i n tim e with the modified Kawaler spin down model 
described in i|2.4.1l The black curve shows the critical rotation 
rate (ui C rit at 550 Myr) above which stars fall in the saturated 
regime. The dashed curves show the median rotation period at 
each age. The associate mass and spectral type are given along 
the top axis. Masses are taken from a 1 Gyr, [Fc/H]m37=+0.045 
YREC isochrone. Convergence occurs later for redder (lower mass) 
stars. 
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analogy to F igure [3] M37's color-period distri- 
bution (Hartman ct al. 2009) evolved forward in time from 0.55 
Gyr (red stars) to 10 G yr (purple crosses) accor ding t o the 
IReiners fc Mohantvl 1 120121) braking law described in H2.4.2I The 
corresponding ages (in Gyr) are indicated to the left of the mass- 
rotation sequence. The black line shows the constant critical ro- 
tation rate (t*; cr i t =8.56xl0 -6 s — 1 ). The dashed lines indicate the 
corresponding medians from the Kawaler wind law (Figure [3j for 
comparison. The associate mass and spectral type are given along 
the top axis. Masses are taken from a 1 Gyr, [Fe/H]M37=+0.045 
YREC isochrone. 
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Gyr. Even in the initial cluster distribution at 550 Myr, 
bluer stars have begun to form a mass-rotation sequence. 
Because of the rotation-dependent rate of angular mo- 
mentum loss (Eq. |4|) , this sequence becomes tighter over 
time and redder stars converge to join it. 

The vertical displacement of the coeval distributions 
reflects the fact that stars spin more slowly with time. 
The reddest stars in the cluster exhibit a large dispersion 
in rotation rate at young ages, but this range also shrinks 
over time. Nonetheless, a star with spectral type M4 and 
a rotation period of 20 days easily falls within the over- 
lapping range of the 0.55, 1, and 2 Gyr distributions, but 
could also be among the most rapidly rotating members 
of even older populations. This inherent ambiguity rep- 
resents a fundamental uncertainty in rotation-based ages 
and will be discussed in depth in £ 13.2. 11 

The median mass-rotation curves at a series of ages 
may be inverted to provide rotation-based ages when the 
scatter about the median is small. However, we note 
that rotation rates still encode age information, even in 
situations when a star's mass and rotation rate does not 
correspond to a unique age. If the form of the distribu- 
tion was known, and the rate of change calibrated, one 
could in principle obtain statistical information about 
ages even when the rotation-age solution is not unique. 
A probability distribution for the age of a star with a 
given mass and rotation period could be constructed by 
weighting these iso-age rotation distributions by the rate 
of star formation and the rate of angular momentum evo- 
lution. Given the uncertainties associated with the rates 
necessary to this calculation, we instead focus on the 
median rotation-mass-age relationship here and use the 
range of rotation to quantify the uncertainty in rotation 
based ages in §3.21 

For comparison, we also explore the the effect of using 
a different angular momentum loss law. We begin with 
the same M37 ro tation distribution at 550 Myr, but this 
time we use the iReiners fc Mohantvl wind law ( H2.4.2p 
to model angular momentum evolution. Figure [4] shows 
how the shape of the M37 distribution changes under the 
IReiners fc Mohantvl wind law compared to the modified 
Kawaler wind law (Figure |3|) . Conspicuously, the slope 
of the mass-rotation relationship appears to flip between 
550 Myr and 1 Gyr. This arises from a combination of 
negligible spin down for low mass stars and an interval 
of rapid spin down for solar-type stars as noted in Figure 
[2l After 1 Gy r, the rate of angula r momentum loss pre- 
dicted by the IReiners fe Mohantvl law slows for stars of 
all masses, maintaining the inverted slope of the mass- 
rotation sequence established earlier. 

The dashed lines in Figure [4] reproduce the median 
mass-rotation sequenc es from the Kawale r law. The 
difference between the IReiners fe Mohanty and Kawaler 
sequences is striking In addition to ha ving slopes of 
a different sign, the IReiners fc Mohantvl law produces 
a smaller torque than the Kawaler law. To wit, the 
Kawaler law spins M < 0.9 M Q stars down to rota- 
tion pCTiods_sk2wer_than they achieve in 10 Gyr with 
the IReiners fc Mohantvl law- in half the time! Because 
of the compressed dy namic range of rotation rates in the 
IReiners fc Mohantvl law compared to Kawaler, the over- 
lap between iso-age populations is more extreme for low 
mass stars. The reddes t stars have suc h weak predicted 
spin down with a Rei ners fc Mohantvl law that the 550 
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Fig. 5. — Mass-period distribution for M37 sample 

IHartman et al.l H2009D ■ The dark gray band marks the smoothed 
25 th and 75 th period percentiles. The light gray band depicts the 
smoothed 10 th and 90 th period percentiles We generally adopt 
the interquartile region as the proxy for width of the rotation 
distribution, although the im pact of using the inner 80% region 
instead is discussed in i|3,2,3l 



Myr and 10 Gyr distributions overlap. In this regime 
where spin down failed to separate iso-age populations, 
the overlap creates ambiguous rotation-based age predic- 
tions covering a wide range. To extract more informa- 
tion, one must rely on statistical ages. 

Now that we have two predictions for how M37's rota- 
tion distribution evolves in time, we need to characterize 
their median trend and range. To quantify the width 
of the rotation distribution, we adopt the interquartile 
range in period. We calculate the 25 th and 75 th per- 
centile in a moving window of 25 stars per mass bin and 
perform a smoothed boxcar average with a width of 50 
stars to flatten out any star-to-star jitter in the calcu- 
lated percentiles. The middle 50% of stars in the ob- 
served M37 cluster mass-rotation distribution lie within 
the interquartile band shown in Figure [5] and we use this 
to characterize the spread between fast and slow rotators 
throughout most of the paper. The median is calculated 
with the same bins and smoothing as for the interquartile 
region. 

For comparison, we also consider the 10 th and 90 th per- 
centile rotation periods because they set much broader 
bounds than the interquartile region, but are not overly 
sensitive to outliers in the rotation distribution. We cal- 
culate the 10 th and 90 th percentile in a moving window 
of 25 stars per mass bin and perform a smoothed boxcar 
average with a width of 80 stars. The middle 80% of 
stars in the observed M37 cluster mass-rotation distri- 
bution fall within the light gray band shown in Figures 
[5) The effect of choosing this alternative metric for the 
width of the rotation distribution is discussed in ^3.2.31 

The percentile ranges in Figure [5] clearly illustrate that 
higher mass stars have converged to a narrow mass- 
rotation sequence by the age of M37. Toward lower 
masses, the distribution rapidly fans out to span over 
a decade in period. This sharp transition may be used 
to define the point of convergence for which more mas- 
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sive stars have forgotten their initial rotation rates and 
converged to a tight sequence. We choose to define the 
point of convergence as the place where the width of the 
distribution is comparable (within a factor of 2) to the 
measured period uncertainties. As we shall see later, sur- 
face differential rotation limits the precision of rotation 
period measurements. In H3.2.21 we adopt a minimum 
uncertainty of 10% to reflect the range of surface rotation 
rates at different latitudes. Thus, we define our criterion 
for convergence to be when the interquartile region is 
smaller than 20% of the median rotation period. This 
also qualitatively reproduces the by-eye location of the 
abrupt widening of the distribution. 

By this metric, a converged sequence is established 
down to (B-V)o =1.08 (M=0.78 M ) at the starting 
age of 0.55 Gyr. For the Kawaler wind law, the mass- 
rotation sequence is converged to (B-V)o =1.17 (0.74 
M Q ) by 1 Gyr, to (B-V) =1.41 (0.62 M Q ) by 2 Gyr, 
to (B-V)o =1-50 (0.52 M Q ) by 5 Gyr, and to (B- 
V)o =1-54 (0.46 M ? ) by 10 Gyr. By the same rule, the 
iReiners fc Mohantvi distributions are converged to (B- 
V) =1.19 (0.73 M Q ) by 1 Gyr, to (B-V) =1-33 (0.67 
M ) by 2 Gyr, to (B-V)o =1.46 (0.58 M©) by 5 Gyr, 
and to (B-V)o =1.54 (0.45 M Q ) by 10 Gyr. It is inter- 
esting that the location of the point of convergence is 
robust to within 0.06 M Q even though these two models 
treat the saturation threshold differently. These points 
define the red boarder of the regime where the distri- 
bution of rotation rates has tightened enough that the 
median rotation-mass-age relationship can provide good 
constraints on age under one of these wind laws. 

3.1.2. Comparison to Published Gyrochronology 
Relationships 

In the previous section, we quantified both the me- 
dian rotation-mass-age relationship and range of rota- 
tion rates for two different angular momentum loss laws. 
Here, we reconcile those trends with the literature. Pub- 
lished gyrochonology relationships ignore rapidly rotat- 
ing stars and focus on describing the shape and time- 
evolution of the sequence of slowly-rotating stars in open 
cluster samples (for further discussion, see i B.2l and ref- 
erences therein). For the purposes of comparison, we 
concentrate on the median trend and postpone further 
discussion of the width of the rotation distribution until 
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Figure |6] shows the median mass-rotation relation- 
ship obtained by evolving the M37 rotation distri- 
bution forward wit h both a modified Kawaler and 
IReiners fc Mohantvl wind law (from <j3.1.1l) . It also 
features the p rediction s of three g yrochronolo gy re- 
lations (IBarnesI l2007t iMamaiek fc Hillenbr and 2008; 
iMeibom et alJ 12009). with the functional form of Eq. []] 
and parameters specified by Table [TJ 

The medians show a steep upturn at low mass that 
is not present in the gyrochronology curves. This under- 
scores a fundamental difference: the median is a measure 
of the entire distribution, while published gyrochronol- 
ogy relationships are calibrated only to the slow rotators. 
A comparison between these two is fair in the converged 
regime because there the entire rotation distribution has 
collapsed to the narrow, slowly-rotating band fit by gy- 
rochronology relationships. Figure [6] may be interpreted 
by assuming everything to the left of the upturn in the 
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Fig. 6. — Median mass -period relationship for a Kawaler (black 
solid) and Rcincrs & Mohanty (black dotted) wind law at fixed 
age ranging from 0.55 Gyr (a) to 5 Gyr (d). Isochrones 
for three gyrochronology relations (Barnes 2007, dash-dot; 
Mamajck & Hillenbrand 2008, short dash; and Mcibom ct al. 2009J 
long dash) are shown in gray for comparison. For reference, the 
Sun is shown at its current equatorial rotation period in all panels. 
The color-mass conversion used YREC isochrones of the respective 
age and [Fe/H]=+0.045. 



median trend falls in the converged regime (see ^3.1.1l for 
a more formal definition of the point of convergence.) 

Although open clusters provide information about ro- 
tation in the young stars, the Sun is the only point used 
to calibrate angular momentum evolution models at ages 
older than 1 Gyr. Given its significance, we first check 
how well these models reproduce the Sun. For a vi- 
sualization, panel (d) comes the closest in age to the 
Sun at 5 Gyr. By definiti on, the modified Kawaler and 
IReiners fc Mohantvl (|2012f ) medians reproduce the Sun at 
the solar age. Two of the three gyrochronology relations 
also match the Sun's equatorial rota tion rate within ±1 
day. The third. IMeibom et a l. (2009), predicts a rotation 
period ~ 20 days for the Sun be cause they did not use it 
as a calibration point. In stead, Mcibom et al.l retained 
the time-dependence from IBarnesI ()2007l ) while changing 
the color-rotation function f(B — V) to fit the open clus- 
ter M35. Nonetheless, the gyrochronology relations best 
agree with our median rotation-mass-age relations near 
solar color. 

For stars redder than the Sun, all three gyrochronol- 
ogy relations spin stars down too quickly, resulting in 
predicted rotation periods that are too long by a factor 
of 1.5 — 2 relative to a modified Kawaler spin down model. 
These published relations display zero-point differences 
relative not only to theory (panels b-d), but also to the 
observed rotation distribution of M37 (panel a). We note 
that the coupling timescale i s shorter than the age of 
M37 (Dcnisscn kov et aT1l2010l ). so this cluster's rotation 
distribution should reflect the true angular momentum 
content of its stars. Gyrochronology relationships may 
predict rotation rates that are too slow because they are 
calibrated to younger clusters with stars whose core and 
envelope have not yet re-coupled. If the surface rota- 
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tion rate masks a reserve of angular momentum in the 
core, gyrochronology relations would propagate forward 
the superficial rotation rate and predict artificially slow 
r otation at late times. 

IMamaiek fc Hillenbrand! (|2008f ) noticed that the 
Barnes (2007) relation predicted rotation periods 50% 
longer than those found in the Hyades. To address this, 
IMamaiek fe Hillenbrand! carefully selected only stars 
that clearly lie on the converged sequence to include in 
their fit to t he ma ss-rotation function f(B — V). As 
a result, the [Mamajek & Hillenbrand isochrones repro- 
duce the shape of the mass- rotation sequence better than 
other gyrochronology prescriptions. Differences between 
the shape of the gyrochronology curves and the median 
of the M37 data in panel (a) arise strictly from problems 
in the calibration of gyrochronology relations in young 
and intermediate age clusters. However, additional time- 
dependent differences can also arise because our angular 
momentum models evolve in time and the gyrochronol- 
ogy relations remain self-similar. Two effects in partic- 
ular can impact the rotation distribution and therefore 
its median: evolutionary changes to stellar structure and 
magnetic saturation. The former is most important for 
the rapidly evolving, relatively high mass stars. The pre- 
scriptions for J also depend on the mass and radius. This 
structural dependen ce is much more signficant for the 
iReiners fc Mohantyl law than for the modified Kawaler 
law. Magnetic saturation is important for both loss laws. 
In the case of the modified Kawaler law, only low mass 
stars rotate faster than the saturation threshold at the 
age of M37 because of the Rossby-scaling (see Figure 
[3] and Eq. [5]). In contrast, the sat uration threshold is 
constant in the IReiners fc Mohantvl law and affects stars 
of all masses because of the way it cuts through M37's 
rotation distribution (Figure 2]). Together these two ef- 
fects can have a mild (modified Kawaler) to significant 
(|Reiners fc Mohantvf ) time-dependent impact on the me- 
dian rotation-mass-age relationship. 
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Tests in Old Field Stars: The Importance of 
Differential Rotation 



3.1.3. 



Older than 1 Gyr, the rotation-mass-age relationships 
described in ^3.1.21 are constrained only by sparse data, 
primarily calibrated on the Sun. We are fortunate to 
possess rotation period and age measurements for two 
wide binaries. These old field stars can test whether the 
solar-calibrated rotation-mass-age relationship holds for 
stars with a mass and age different than the Sun. 

Because stars lack a uniform surface rotation rate, any 
rotation period measurement of a star will be drawn from 
some distribution of surface rotation periods. If you mea- 
sure a rotation period using a near-polar starspot, it will 
yield a different rotation period relative to an equato- 
rial s tarspot because of latitudinal di fferential rotation 
(e.g. iBarnes et all [20051: iRemersI [20061 1 The Sun ex- 
hibits differential rotation, with a rotation period of 25.4 
days at the equat o r and 36 da ys around its poles (e.g. 
iSchou et all 119981 : iBeckl 120001 and references therein). 
Disk integrated Ca II K-line measurements yield a range 
of rotation rates from 24.5 (lower than the true equato- 
rial rotation rate due to me asurement error) to 2 8.5 days, 
with a mean of 26.09 days (|Donahue et al.|[l996l) . Below, 
we investigate the data available for the two binary star 
systems: 61 Cyg and a Cen. 



iDonahue et al.l (|1996l ) present the measurement of a 
range of rotations for 36 stars with five or more seasons of 
detectable periodic modulation in the Ca H and K emis- 
sion from the Mount Wilson survey. Among these, the 
K5V+K7V binary 61 Cyg (HD201091 and HD201092) is 
of particu lar interest beca use it is physically well char- 
acterized. IDonahue et al.l report mean rotation periods 
{P) A = 35.37 days and (P) B = 37.84 days with pe- 
riods ranging from 26.82-45.25 days for 61 Cyg A and 
31.78-46.57 days for B. The surface metallicity is spec- 
troscopically measured as [Fe/HU = —0-20 ±0.10 and 
[Fe/H] B = -0.27±0.19 (|Luck fc Heitedl2"005h . The high- 
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accuracy parallax information from 
1995) and the Hipparcos catalog (jPerrvman fc ESA 



1997) tightly constrains the luminosity to La — 0.153 ± 
0.010 L Q and L B = 0.085 ± 0.007 L . Combining 
the parallax w i th int erferomctric radius measurements, 
IKervella et all (2008) found R^ = 0.665 ± 0.005 R Q and 
R_b = 0.595 ± 0.008 R . T wo measurements o f the 
mass exist: 0.67 a nd 0.59 M^ dWalker et al.l ll995l) and 
0.74 and 0.46 M Q (IGorshanov et al.1120061) for 61 Cyg A 
and B, respectively. IKervella et al.1 (J2008I ) used observa- 
tional c onstraints to construct CESAM2k evolutionary 
models (|Morellll997h and derive masses M^ =0.690 and 
M B =0.605 M and an age of 6 ± 1 Gyr. 

We adopt the IKervella et al.l masses and re-derive the 
age with YREC models. We evolved models using the 
input physics described in van Sa ders fc Pinsonneaultl 
(2012), adopting a solar calibrated mixing length of 
a = 1.926 and Y = 0.265. The initial metallicity was 
chosen to yield a surface [Fe/H] of —0.20 at an age of 6 
Gyr. We matched the radius within the errors at ages of 
5.7 ± 1.2 and 7.5 ± 2.5 Gyr for 61 Cyg A and B respec- 
tively. This yields an age of 6±1.4 Gyr for the system as a 
whole, to be compared wi th the age of 6 ± 1 G yr derived 
using the CESAM code (jKervella et al.ll2008T ) with the 
same masses and helium, but different choices of metal- 
licity and mixing length. However, we believe that this 
underestimates the age error, as an additional 1 Gyr sys- 
tematic results from varying the metallicity with the ob- 
servational error of 0.1 dex. We therefore adopt 6.0± 1.7 
Gyr as our age for the system. 

a Cen 

The proximity of a Cen allows the precise determina- 
tion of the fundamental parameters for its component 
G2V+ K1V stars. We adop t ed th e same physical param- 
eters as lEggenberger et al.l ([2 004). These include precise 
parallaxes (|Soderhichn 1999), masses (|Pourbaix et al.l 
2002) , radii (LK crvclla et al.l I2003D. and temperature s 
(jNeuforge-Verheecke fc Magainlll997l : lMorel et"aTll2000( ). 
Detections of p-mode oscillations in both a Cen A and 
B provide seismic constraints on the central helium 
abundance and age of th e star (|Bouchv fc Carrier] 120021 : 
iCarrier fc Bourbanll2003D . Combi ning all of these non- 
seismi c and seismic constraints, Mig lio fc Montalbanl 
(2005) found an age of 6.8 ± 0.5 Gyr for a Cen. As 
described above, we created a YREC model for stars of 
mass M A =1.105, M B =0. 934 and Z = 0.0302, Y = 0.275 
(|Eggenberger et al.|[2004D . 

Because a Cen is not observable from the northern 
hemisphere, long-term monitoring of rotation rates from 
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Fig. 7. — Location of the Sun and wide binaries a Cen and 61 
Cyg in the effective temperature-rotation period plane. The Mount 
Wilson survey measured the range of observed rotation periods over 
many seasons due to surface differential rotation (dashed), mean 
rotation rate (point), and erro r in the mean (solid) for the Sun and 
61 Cyg (Donahue ct al. 1996). Because a Cen lacks such detailed 
rotation measurements, we show the two measurements available 
for a Cen A and assume that a Cen B displays a range in rotation 
rates similar to 61 Cyg A, which has a similar rotation period. The 
gray contours represent the median temperature-period relation for 
stars with age 1-10 Gyr, assuming a modified Kawaler wind law 
and YREC isochrones. The hatched regions illustrate the region 
consistent with the isochrone ages for the two binaries. 

th e Mount Wilson survey are no t available. We refer 
to lMamaiek fe Hillenbra nd (2008), which cites two mea- 
surements of rotation period for a Cen A: 22 ± 3 days 
(Gui nan 2008, private c ommunication) and 28.8 ± 2.5 
days (|Hallam et al.lll99l"l ). The rot ation period of a Cen 
B is measured as 36.9 ± 0.5 days (|Jav et al.lll997[ ) from 
the 1995 IUE data set (a more recent analysis of the 
same d ata yields a nearly identical result 36.2 ±1.4 days, 
DcWarf ct al. 2010). Because a Cen B has only one inde- 
pendent measurement of its rotation period, we assume 
that its surface differential rotation produces a range of 
rotation rates as large as 61 Cyg A. This is a reasonable 
assumption since a Cen B and 61 Cyg A poss ess simi- 
lar rotation periods and IDonahue et al.l (1996) found a 
strong correlation between mean rotation period and the 
range of observed rotation rates. 

Surface Differential Rotation in Binaries 

Although individual rotation periods may be measured 
precisely, they are drawn from a range of rotation periods 
corresponding to the latitude of the active region on the 
stellar surface. The range of rotation periods observed 
for the Sun and 61 Cyg and inferred for a Cen are shown 
in Figure [7] These ranges are compared with the median 
rotation curves evolved forward in time from the M37 
rotation distribution assuming a modified Kawaler wind 
law. 

Because binaries form coevally, the rotation rates of in- 
dividual components should lie along the same rotational 
isochrone. Furthermore, the shaded regions indicate the 
range of rotation periods consistent with the seismic age 
of a Cen A and isochrone age of 61 Cyg. The measured 
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Seismic or isochrone Age (Gyr) 

Fig. 8. — Comparison of rotation-based age with isochrone age 
for the Sun and wide binaries a Cen and 61 Cyg. Rotation based 
ages calculated according to the median temperature-rotation re- 
lationship. As in Figure the observed range of rotation rates 
are measured for the Sun and 61 Cyg. Since a Cen B has a single 
rotation period measurement, its range is assumed to be of compa- 
rable size to 61 Cyg A. Two individuals measurements of rotation 
period are shown for a Cen A, which together can be taken to set 
a minimum on its observed range of rotation rates. For clarity, 
the central age values of 61 Cyg B and a Cen B have been offset 
slightly from their true isochrone age. 



rotation periods for 61 Cyg A cover a broad range of 
ages from approximately 3-9 Gyr. The hatched region 
corresponds to the range of rotation periods expected for 
a star with age 6.0 ± 1.7 Gyr. Both components of 61 
Cyg are consistent with the isochrone age and with each 
other. 

For a Cen, the tightly constrained seismic age is con- 
sistent with t he B component and the longer period 
measurement (iHallam et alj[l99ll) of A's rotation period. 
The short period measurement of a Cen A rotation com- 
pared to the isochrone age demonstrates that any single 
measurement can have a ~ 2 Gyr error in its rotation- 
based age due to surface differential rotation. 

To better quantify the range in rotation-based ages due 
to a range in observed rotation rates, we model the angu- 
lar momentum evolution of these stars using a modified 
Kawaler wind law and the customized stellar evolution- 
ary tracks described above. The initial rotation period is 
set to M37's median rotation period at the correspond- 
ing temperature and is evolved forward from 550 Myr. 
We use temperature as the coordinate to compare the 
near-solar metallicity M37 rotation distribution with the 
binaries, which have different metallicities and therefore 
different mass-color relationships, because the convective 
properties that govern angular momentum evolution are 
primarily set by T e fi- We interpolate on this track to de- 
termine the rotation-based age corresponding to range of 
observed rotation rates. The results are compared with 
the seismic age of a Cen or the isochrone age of 61 Cyg 
as shown in Figure [8] 

Because these ages were computed with customized 
tracks, they will differ from the ages read off the near 
solar-metallicity M37 contours shown in Figure [Jj In 
particular, this drives rotation-based ages for the slightly 
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metal-poor and lower mass 61 Cyg into worse (la- 2a) 
agreement with its isochrone age, which hints at a possi- 
ble metallicity dependence to the angular momentum loss 
law. Because the only point of calibration for loss laws 
at ages > 1 Gyr is the Sun, it is conceivable that stars of 
a different metallicity or mass would behave differently 
beyond the regime tested by clusters. The Kepler mis- 
sion will create a sample of field stars with asteroseismic 
ages capable of calibrating the metallicity dependence of 
the rotation-mass-age relationship. 

We also perform an id entical calculation using the 
iReiners fe Mohantvl (120121) wind la w and the customized 



tracks. The Rei ners fe Mohantvl models do not spin 
down to periods longer than 30 days for a Cen B and 
61 Cyg A and B within 10 Gyr. This model then fails to 
provide a rea sonable rotation-based age for these binary 
components. iReiners fe Mohantvl (|2012l ) note a related 
problem for field M stars, which are below the mass range 
considered here. However, even for K dwarfs slowly ro- 
tating stars exist with well measured ages younger than 
10 Gyr. This indicat es that the late time behavior of the 
IReiners fe Mohantvl law as published does not appropri- 
ately capture true stellar behavior even for higher mass 
dwarf s. Further discussion of the IReiners fe Mohantvl 
(|2012f) la w appears in $£21)1 

If the IReiners fe Mohantvl rotation-based ages are 
too old, some published gyrochronology relationships 
produce ages that are too young. We compare 
the two solar-calibrated published gyrochronology re- 
latio nships against these bin a ries a nd determine that 
the iMamaiek fe Hillenbrand! ([20081 relation is con- 
sistent while the Barnesl (120071 ) formalism is not. 
lMamaiek~ fe Hillen brand! (J2008D already performed this 
comparison, finding an age (in years) of log tA — 9.57 and 
log is = 9.53 of 61 Cyg's components using their mean 
rotation period, which average to (logi6iCyg) = 9.55. Al- 
though this mean age (3.52 Gyr) is too young compared 
with the isochrone age of 61 Cyg (6.0 ± 1.7 Gyr), includ- 
ing the full range of rotation observed in 61 Cyg broad- 
ens the range of rotation-based ages to 2.25 < tA < 5.66 
Gyr and 2.48 < in < 4 .87 Gyr. Similarly, for a Cen, 
iMamajek fc Hillenbrand! report (logt Q ccn) = 9.70 (5.01 
Gyr), which is inconsistent with the seismic age (6.8±0.5 
Gyr) at the 3.6a level if the gyroage is treated as exact. 
With the inferred range of rot ation periods, the tension 
between the seismic age and the IMamajek fc Hillenbrand 
age ranges (3.64 < t A < 5.86 Gyr and 4.77 < t B < 
5.47 Gyr) decreases to the ~ 2a level. Thus, the 
mean IMamaiek fc H illenbrand trend is consistent with 
these binaries, but their estimated errors are too small. 
IMamaiek" fc Hillen brand! take the difference in predicted 
rotation-based ages of these coeval binary components 
as the error. They find an uncertainty of ±11% in gy- 
rochronol ogy pred i ction s, similar to the total error of 
15% from Barnes (2007). These error estimates ignore 
systematic sources of uncertainty in rotation-based ages, 
including surface differential rotation, which we will dis- 
cuss in §3.21 

We next consider the Barnes (2007) relationship. For 
conciseness, we write the age corresponding to the mean 
and range due to surface differe ntial rot ation in a con- 
densed format. For 61 Cyg, the iBarnes! relation yields 
ages of t A = 2.52±i;q4 Gyr and t B = 2.02to;58 Gyr. For 



a Cen, we find t A = 4.67±"i Gyr and t B = 4.29±g;|f 
Gyr. Even with the generous inferred range of surface 
rotation due to differential rotation, the IBarnes! relation 
predicts ages that are systematically too young, espe- 
cially in the case of a Cen B. This large discrepancy with 
old field binaries, coupled with its deviati on from M37's 
median trend (Figure[f)]), indicates that the Barnes (2007) 
relationship does not reproduce observations. Better ro- 
tation measurements for the a Cen system to determine 
a true, rather than inferred, range of differential rotation 
are important for testing rotation-mass-age relationships. 
These examples underscore the fact that differential ro- 
tation imposes an uncertainty in rotation-based ages and 
we will return to this point in §3.2.21 

In addition to these binary systems, photometric mon- 
itoring programs have reported rotat ion periods for large 
samples of FGK (CoRoT lAffer et all I20TI ) and M 
(MEarth. Hrwin et al.ll2011| ) field stars. Although these 
represent a heterogeneous stellar population that will in- 
clude synchronized binaries as well as both dwarfs and 
subgiants, the rotation distribution gleaned by these 
studies still provides interesting constraints. These field 
star studies report rotation periods as long as 80-100 
days for M-F stars. The existence of such slowly spinning 
stars presents a challenge to current rotation-mass-age 
relationships, if they are main sequence stars. Because 
subgiants are expected to have longer rotation periods 
than dwarfs, we recommend a cautious interpretation of 
field star observations that do not distinguish evolution- 
ary state. The three gyrochronology relationships tested 
here all predict that a 100 day period corresponds to ages 
older than 10 Gyr within the color regime B-V< 1.5. 
IReiners fc Mohantvl (|2012f ) recognize that their wind law 
is m uch too weak to exp lain the slowly rotating M stars 
from Irwi n et all (J2011). To account for this problem, 
they change their saturation threshold from 8.5 days to 
P crit = 40 days below M < 0.3 M . They justify this 
break in the saturation behavior with an argum ent that 
the effecti ve tcz rises sharply at low masses (e.g.lSt epien 
Il994t iKim fc Demarque! 119961: iKiraea fc S tcpicn 2007). 
With this adjustment, the Rei ners fc Mohantvl law re- 
produces the ~ 100 day rotation periods for M stars, but 
retains rapid rotaion fo r more massive stars . Even with 
the adjusted P C rit> the IReiners fc Mohanty law cannot 
explain the recent detec tion of equally lo ng rotation pe- 
riods in K and G dwarfs (jAffer et al.ll2012l ). The modified 
Kawaler wind law also fails to predict 100 day periods. 
At 10 Gyr, our models predict a median rotation period 
of ~ 40 days for G dwarfs, dropping to ~ 70 days for 
M dwarfs. These long-period field stars will present a 
serious challenge to our understanding of angular mo- 
mentum evolution if the period measurements prove to 
not be aliased and if further examination places these 
stars on the main sequence. 

3.2. Uncertainties in Rotation-Based Ages 

Until stars lose the memory of their initial conditions, 
we cannot distinguish between an inactive young star 
and an active older star at the same rotation period. 
This implies that there exists a minimum error in the 
mean rotation-mass-age relationship due to the intrin- 
sic dispersion in stellar rotation rates. We quantify this 
minimum age uncertainty and evaluate in what mass and 
period regimes the rotation period becomes a good esti- 
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mator of age for field stars ( H3.2.1[) . 

Furthermore, rotation period measurements suffer 
from uncertainties. While instrumental and observa- 
tional errors may be removed by longer baseline monitor- 
ing or higher precision photometry, the fact that stars do 
not have a unique surface rotation period presents an ir- 
reducible uncertainty (s ee ij4.1 | for a possib le exception). 
Seismic inversion (e.g. iSchou et al.l 119981) and observa- 
tions of s unspots at d ifferent latitudes on the Sun's sur- 
face (e.g. [Lustig 1983) agree that the solar rotation rate 
is significantly faster at the equator than near the poles. 
Because we cannot know if we are measuring equatorial 
rotation rates for unresolved stars, this sets bounds on 
the precision of rotation period, and therefore age mea- 
surements using this technique ( £|3.2.2[) . 

We also consider how different modeling choices influ- 
ence our derived uncertainties. Section [3 . 2 . 31 investigates 
the impact of choosing a broader metric to describe the 
spread in the rotation distribution. The uncertainty in 
the age of M37, which sets our initial conditions, intro- 
duces a systematic error in the age scale ( £13.2.41) . Due 
to uncertainties in color-temperature relationships, dif- 
ferent stellar isochrones yield divergent mass scales for 
a cluster main sequence, particularly among red stars. 
Because the angular momentum evolution depends on 
the internal stellar structure, the adopted mass s cale af - 
fects the evolution of the rotation distribution ( E| 3 . 2 . 5 1) . 
Finally, we examine the impact on rotation-based age un- 
certainties of using the two different angular momentum 
loss laws described in 



3.2.1. Uncertainty Due to Initial Conditions 

For our base case, we adopt a modified iKawa lcr (1988) 
angular momentum loss model. We distill the M37 dis- 
tribution into a mass-dependent median rotation period 
and use the interquartile region to measure of the period 
uncertainty due to the spread of initial rotation rates. 
We evolve the smoothed 25 th and 75 th percentile curves 
(defined in £13.1. 1|) forward in time until we reach ei- 
ther the age of the Galactic disk, which we take to be 
t = 10 Gyr, or when a model corresponding to that mass 
leaves the main sequence. Table [5] contains the 10 th , 
25 th , 50 th , 75 th , and 90 th percentile curves as a func- 
tion of time, where (B-V)o and T e gf are computed at 
fixed mass using a YREC isochrone of the appropriate 
age and [Fe/H]=+0.045. Interpolation on this grid al- 
lows the determination of rotation-based ages and the 
associated uncertainty. 

The top panel of Figure |H] shows an example of how the 
interquartile range changes over time at fixed mass. A 
horizontal cut across these curves gives the 25 th and 75 th 
period percentiles for a bin centered at M = 0.60 M at 
fixed age. In this calculation, the predicted range of ro- 
tation periods converges to a unique value at fixed mass 
and age. For field stars, it is useful to invert the problem 
and ask what range of ages is expected for a measured 
rotation period. This corresponds to taking a vertical cut 
through Figure [SJ which gives the range in age predicted 
for a fixed rotation period. A unique rotation-mass-age 
relationship will predict a spurious range of ages for co- 
eval stars born with different initial rotation rates be- 
cause stars that are initially rotating faster take longer 
to reach a given period than initially slow rotators do. 
We take the difference between ages inferred from the 
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Fig. 9. — Top: Interquartile range of periods from M37 evolved 
forward using a solar calibrated Kawaler law for a bin of 50 stars 
centered on 0.60 Mq (black lines). Vertical red lines cut the rota- 
tion profile at fixed period showing range of possible ages. Arrows 
indicate that gyrochonology can place either an upper limit on the 
age of a star with this period, but that it may be younger than 550 
Myr; or, a lower limit, bounded by the age of the Galaxy. Bottom: 
Width of interquartile band as a function of period. The quantity 
A log(Age) gives the uncertainty in rotation-based age. 



top and bottom of our rotation distribution (A log Age) 
as the uncertainty in the age estimate from gyrochronol- 
ogy (Figure \§§. This ambiguity in age is a fundamen- 
tal limitation of this technique; it takes a finite time for 
stars to forget their initial conditions. The lower panel 
shows the size of the age uncertainty for different mea- 
sured rotation periods at fixed mass. The uncertainty in 
age declines toward longer periods because, if stars have 
time to spin down to long rotation periods, the rotation 
distribution will have had more time to converge. 

There exists a region where gyrochronology can only 
provide an upper limit on the age of a star. We define 
this to be the case if the measured period lies within the 
interquartile range of M37 at 550 Myr. We lack further 
information due to our choice of initial conditions, and 
we discuss rotation-age distributions in young stars in a 
companion paper (Paper II). Only ~ 5% of field stars 
should fall in the regime where only an upper limit is 
relevant, assuming the Galactic disk to be 10 Gyr old 
with a constant star formation rate. On the other end, 
at 10 Gyr, the period distribution still has some width 
and the slowly rotating quartile period sets a lower bound 
on the age and the upper bound is the age of the Galaxy. 

Expanding our focus, we determine how the uncer- 
tainty in age depends not only on rotation period, but 
also on mass. In this ideal case, where we have assumed 
both perfect mass and rotation measurements, and that 
the theory fully encapsulates the physics of stellar an- 
gular momentum loss, the contours on the top panel of 
Figure [10] give the best possible constraints on age for 
each combination of mass and period, where lighter grays 
correspond to smaller A log Age. The uncertainty in age 
increases towards lower mass because convergence takes 
longer for these stars. Gyrochronology can define an age 
range, rather than just place upper limits on the age, in 
the shaded region. In some windows of mass and period 
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TABLE 2 
Median Gyrochronology with Uncertainties for a Modified Kawaler 

Wind Law 



Age 


Mass 


(B-V)o 


Toff 


PlO 


P25 


P50 


P75 


P90 


(Gyr) 


(M Q ) 


(mag) 


(K) 


(days) 


(days) 


(days) 


(days) 


(days) 


0.55 


0.45 


1.54 


3643 


0.58 


0.77 


2.09 


9.43 


14.42 


0.58 


0.45 


1.54 


3644 


0.59 


0.79 


2.14 


9.65 


14.77 


0.63 


0.45 


1.54 


3643 


0.62 


0.83 


2.25 


10.15 


15.53 


0.69 


0.45 


1.54 


3642 


0.66 


0.87 


2.38 


10.71 


16.39 


0.76 


0.45 


1.54 


3642 


0.69 


0.92 


2.52 


11.35 


17.36 



Note. — Table 2 is published in its entirety in the electronic edition of this 
journal. A portion is shown here for guidance regarding its form and content. 
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Fig. 10. — Uncertainty in age (A log(Age/Gyr)) as a function of 
measured period and mass. The shaded region shows where gy- 
rochronology is valid for our models. Uncertainties in age are cal- 
culated in the same manner as for Fig. [9] Lighter grays correspond 
to better constraints on the rotation age. Dotted lines give the me- 
dian rotation period as a function of mass when the M37 distribu- 
tion is evolved to 0.55, 1, 2, 5, and 10 Gyr using a solar calibrated 
Kawaler wind law. The mass axis is converted to temperature 
above assuming a 1 Gyr, [Fe/H]M37 = +0.045 YREC isochrone. 
Both the panels assume that the solar calibrated Kawaler wind 
law perfectly captures the physics of angular momentum loss in 
cool stars and that masses can be measured to arbitrary precision. 
The upper panel assumes that rotation periods can also be per- 
fectly measured, while the lower panel imposes a minimum 10% 
uncertainty in the period measurement. 



space, constraints tighter than Alog(Age) ~ 0.05 dex 
are possible, assuming perfect period measurements. As 
before, the region where gyrochonology can only provide 
an upper limit on the age arises from the initial spread 
in rotation rates. The lower bound occurs because we 
halted the simulation after 10 Gyr. The upward slope of 
the lower boundary at high mass represents these stars 
evolving off the main sequence. We offer the following 
cautions when interpreting this figure. First, because we 
use the interquartile region as our measure of the width 
of the rotation distribution, there will be stars that lie 
outside this range. Second, the specific values of the pe- 
riod and age uncertainty depend on the validity of the 
underlying angular momentum evolution model. 

3.2.2. Impact of Latitudinal Differential Rotation and 
Period Errors 

Unfortunately, achieving these theoretical limits may 
not be feasible with real stars. In ^3. 1.31 we demon- 
strated that old field stars have a significant range in 
measured surface rotation rates which reflects latitudinal 
surface differential rotation. This provides an additional 
fundamental source of uncertainty in rotation-based ages: 
where on the surface is the rotation period being mea- 
sured? 

If the Sun's differential rotation is characteristic of so- 
lar type stars, then the latitude of a star spot could signif- 
icantly impact the measured rotation period. Although 
the range of rotation peri ods may be as large a s 30% of 
the mean rotation period ()Donahue et al.|[l996l) , we con- 
sider a more conservative case where latitudinal differen- 
tial rotation only sets a 10% range in rotation period. To 
account for this, we therefore broaden our interquartile 
range by taking 



A log P = max 



, -P75 cr P 
g P 25 'Plnl0 



(7) 



where crp/P = 10% is the minimum set by differential ro- 
tation. When widening the distribution beyond the ratio 
of quartile periods, we centered AlogP on the geometric 
mean of the quartile rotation periods. We then repeat 
the process described above to obtain the plot shown in 
the bottom panel of Figure [TO] Comparison between the 
upper and lower panels shows that, even given a perfect 
angular momentum loss theory, latitudinal differential 
rotation at this level sets a floor on the age precision 
from gyrochronology estimates instead of allowing the 
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range of rotation rates to converge arbitrarily. 

The impact of surface differential rotation is most se- 
vere in more massive and older stars. For example, if 
a 0.5M Q star has a rotation period of 11.5 days, gy- 
rochronology can only constrain its age to a factor of 4.7 
regardless of whether there is a minimum precision in pe- 
riod. Meanwhile, the age of a 0.8 M star with a period 
of 40 days can be tightly constrained to 2% with perfect 
period measurements, but only to 20% when we include 
a 10% range due to differential rotation. For the Sun's 
mass and equatorial rotation period, its age can be con- 
strained to 2% with perfect period measurements. With 
this correction for differential rotation, the Sun's age is 
predicted to lie between 4.07 and 4.85 Gyr. Generally, 
imposing a 10% range in rotation periods translates to a 
20% uncertainty in age. 

The dotted curves in Figure [10] delineate the median 
rotation period at a selected set of ages. At 550 Myr 
and 10 Gyr, rotation only sets a bound on the stellar 
age. At intermediate ages, the median rotation period 
may be mapped onto a typical uncertainty for a star of 
a given age. Although the logarithmic scale from Fig- 
ure [10] is useful for visualizing the fractional error, the 
absolute age error also provides a helpful comparison. 
The solid lines in Figure [TTJ trace the absolute age un- 
certainties for stars with the median rotation period at 
select ages. Given arbitrarily accurate rotation period 
measurements, stars more massive than 0.75 M always 
have an uncertainty in age smaller than 270 Myr. Adding 
in a 10% uncertainty in rotation period has the effect of 
capping the age precision at around 0.38, 1, or 1.5 Gyr 
for a 2, 5, or 8 Gyr population, respectively. The abso- 
lute age uncertainties are smallest for younger stars, but 
become increasely larger for older stars. Nonetheless, 
the absolute age uncertainty is always below 2 Gyr for 
stars with M > 0.52 M Q . The wide binaries 61 Cyg and 
a Cen corroborate this expected level of uncertainty in 
rotation-based ages due to latitudinal surface differential 
rotation ( £13.1.31) . 

3.2.3. Impact of Choice for Metric for Distribution Width 

In i )3.1.1[ we describe the width of the rotation dis- 
tribution by the interquartile region. One might argue 
that this is too restrictive and does not capture the full 
range of the distribution. We therefore investigate the 
impact of choosing a more generous metric, the range of 
rotation rates between the 10 th and 90 th percentile, to 
characterize the rotation distribution. 

We compare the absolute age uncertainties in Figure 
fTTl The largest differences occur for low mass stars, 
where the spread in rotation rates is maintained until late 
times. The metric chosen to describe the rotation width 
does not change the functional form. Even if arbitrarily 
accurate period measurements are possible, differences in 
the age uncertainty are only a few hundred Myr for stars 
on the tightly converged mass-rotation sequence. From 
the bottom panel, we see these differences due to the 
choice of width metric are smaller than the age uncer- 
tainty due to differential rotation for stars with M > 0.8 
M© at 2 Gyr and to even lower masses for older stars. 
We therefore consider the interquartile region a fair char- 
acterization of the rotation-based age uncertainty. 

3.2.4. Impact of the Age Scale 
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Fig. 11. — Absolute age uncertainty as a function of mass for 
a star with the median rotation period at age 2, 5, and 8 Gyr 
(light to dark). Solid lines correspond to using the inner 50% of 
stars to measure the distribution width while dashed lines indicate 
the inner 80% of stars. The top panel assumes rotation periods are 
measured with arbitrary precision, while the bottom panel includes 
a 10% uncertainty in period measurements. 

Uncertainties in the absolute age assigned to M37 af- 
fect all of our derived rotation-based ages. In §2.1[ we 
initialize all models t o the age of M37 determined from 
lHartman et al.l (12008 1 ). By fitting is ochrones to the clus- 
ter main sequence, lHartman et al.l measured an age of 
550 ± 30 Myr. To see how this uncertainty in the initial 
conditions maps onto our age predictions, consider the 
illustrative case of Skumanich spin down. In that sce- 
nario, the rotation rate scales asw~ £~ a5 . This implies 
that: 



£M37 



^M37 



-i 2 



W* 



(8) 



Thus, the age of M37 (£m37) enters linearly into the 
calculation of a star's age (£*■). An analogous calcula- 
tion may be perfor med for both the Kawaler and the 
iReiners fc Mohantvl wind la ws. As such, the ±5 % un- 
certainty in the age of M37 (H artman et al.1 12008) prop- 
agates to a ±5% uncertainty in the ages predicted by 
these models. This uncertainty in the scale factor repre- 
sents a systematic uncertainty for any semi-empirically 
calibrated model. The random error associated with the 
age of any particular cluster may be reduced by using 
many different clusters. However, the systematic uncer- 
tainty in the overall open cluster age scale will remain. 

3.2.5. Impact of the Mass Scale 

The wind law describes the angular momentum loss 
rate in terms of mass, rather than an observable quan- 
tity like color. Stellar interior models predict divergent 
mass-color relationships, particularly at the cooler end. 
As described in §2.1[ we convert measured (V-Ic)o to 
mass using both YREC and Padova isochrones. Figure 
[U shows that both isochrones agree well for (V-Ic) <l-2, 
but diverge for redder stars. For (V-Ic)>l-2, Padova 
isochrones map stars of a particular color onto a lower 
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Fig. 12. — Absolute age uncertainty as a function of mass for 
a star with the median rotation period at age 2, 5, and 8 Gyr 
(light to dark) calculated with YREC (solid) and Padova (dashed) 
isochrones. The top panel assumes rotation periods are measured 
with arbitrary precision, while the bottom panel includes a 10% 
uncertainty in period measurements. 



mass than YREC isochrones do. Since both isochronc 
systems are widely used for color-temperature conver- 
sions, we investigate how this affects the age uncertain- 
ties. 

We use the Padova isochrones to define the mass scale 
for the M37 initial mass-rotation distribution and repeat 
the calculations described in i )3.2.1l This will have two 
effects on the rotation-mass-age distribution. First, using 
the Padova isochrones should result in a shift of red stars 
to lower masses. Second, the absolute mass scale affects 
the rate of angular momentum loss because it changes the 
star's interior structure, affecting the moment of inertia 
and leading to a different spin down history. 

The absolute age uncertainty derived from the Padova 
isochrones is compared with the YREC results in Figure 
[121 Above 0.77 M , the Padova and YREC isochrones 
produce identical constraints. Below that, the Padova 
isochrones appear to be a stretched version of the YREC 
isochrones. Because the Padova isochrones translate a 
red color to a lower mass, the Padova curves extend to 
lower masses than the YREC ones do. As a result of the 
stretched mass scale, Padova isochrones predict tighter 
age constraints than YREC does for a given mass and 
period. 

Although the exact value of the age uncertainty de- 
pends on the isochrones used to determine the mass scale, 
the difference is confined to the low mass regime and, 
even there, the Padova contours retain the qualitative 
characteristics of the YREC contours. While these sim- 
ilarities are reassuring, it is important to carefully con- 
sider how to map observed colors onto mass. In Figure 
[H it is apparent that the YREC isochrones better repro- 
duce the shape of the cluster main sequence. For this 
reason, we favor the YREC isochrones for this study. 

3.2.6. Impact of the Choice of Angular Momentum Loss 

Law 
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Fig. 13. — Uncertainty in age (Alog(Age/Gyr)) as a function 
of measured period and mass, where the initial mass distribution 
was calculated using YREC isochrones and evolved forward in time 
with the Rcincrs & Mohanty (2012) braking law. The upper panel 
assumes that rotation periods can also be perfectly measured, while 
the lower panel imposes a mi nimum 10% uncertainty in the period 
measurement. See Figure ITul for an explanation of other markings. 



Thus far, we have used the widely-adopted IKawalerl 
(1988) prescription for angular momentum loss, mod- 
ified to include a Rossby-dependent saturation thresh- 
old as described in ^2.4.11 While this simple model re- 
produces observations well in the open cluster domain; 
there is limited data for older field stars. Alterna- 
tive formulations sh ould therefore be explored, and the 
IReiners fc Mohantvl prescription is an intriguing choice. 

Following the pro cedure described in < j3 2.11 with the 
IReiners fc Mohantvl wind law produces the fractional 
age uncertainty con tours displayed in Figure 1131 The 
IReiner s & Mohantv braking law introduces some inter- 
esting features that merit closer examination. Because 
the Rciners law enforces rapid rotation among low mass 
stars until late times, rotation provides only a weak con- 
straint on age. For stars with M> 0.6 M Q , rotation- 
based ages cannot be pinned down to better than 50%, 
regardless of the precision of rotat ion period measure- 
ments. Among higher mass stars, IReiners k, Mohantvl 
predict rapid angular momentum evolution until ~ 700 
Myr. If a star is caught with a rotation period in this 
range, its age may be determined to a precision of about 
7% (or approximately 30-50 Myr). Imposing a 10% pe- 
riod measurement error (Figure 1 131 bottom panel) pro- 



How Good of a Clock is Rotation? 



17 




Mass (M Q ) 



Fig. 14. — Absolute age uncertainty as a function of mass for a 
star with the median rotation period at age 2, 5, and 8 Gyr (l ight to 
dark gray) calculated with Kawaler (solid) and Rcincrs & Mohanty 
(dashed) braking law. The top panel assumes rotation periods are 
measured with arbitrary precision, while the bottom panel includes 
a 10% uncertainty in period measurements. 



duces a sweet spot for stars with rotation periods around 
10 days. 

Figure Q3] shows the absolute a ge uncertainty pre- 
dictcd by both the Kawaler and iReiners fc Mohantvl 
braking laws. When measuremen t uncertainty is in- 
troduced in the bottom panel, the IReiners fc Mohantvl 
law predicts de graded precision in old stars because the 
IReiners fc Mohantvl law has a weaker rate of angular mo- 
mentum loss than the Kawaler law does at ages > 1 Gyr 
(see Figured]). Both exhibit stable rotation-based age 
uncertainties with mass down to a critical point where 
the age uncertainty abruptly becomes quite large, cor- 
responding to the point of convergence. This effective 
boundary to the rotation-based age technique is dramati- 
cally shifted 0.15 M Q higher with the Rcincrs & Mohanty 
law. This shift is du e to the mass- independent way that 
IReiners &: Mohantvl treat magnetic saturation. Compar- 
ison of Figures [3] and [4] shows that IReiners fc Mohantvl 
set the saturation threshold lower than the Rossby scaled 
saturation threshold for (B-V) < 1.4 (M> 0.57 M Q ). 
This keeps stars in the saturated regime longer than with 
a Rossby scaled saturation threshold. Because a Rossby- 
scaled saturation th reshol d is supported by observations 
of open clusters (see l2.4.Ij) . we consider the Kawaler ver- 
sion to be a more realistic boundary to the rotation-based 
age technique. 

Because the rotation-based age constraints differ 
strongly between these two angular momentum loss laws, 
properly testing the rate of angular momentum loss is 
important. For now, we favor the modified Kawaler 
wind law because it properly reproduces Skumanich spin 
down at late time for solar type stars and it incorpo- 
rates a Rossby scaled wind law, both of which are sup- 
ported by open clus ter observations. Additionally, the 
IReiners fc Mohantvl law fails to reproduce the open clus- 
ter distribution, as seen in their Figure 4. They pre- 
dict rotation periods that are too long compared to the 



cluster M50 for M< 0.8 M Q . They argue that this is 
due to core-envelope decoupling, which is not included 
in their modeling. This can explain the presence of slow 
rotators over the mass range where core-envelope decou- 
pling is expected to be important, primarily for solar- 
like stars with shallow surface convection zones. How- 
ever, low mass stars (below 0.5-0.6 M ) are well de- 
scribed empirically with solid body spin down models 
(Tink er~et al.l 12002). The striking differences between 
their models and cluster data for cool dwarfs therefore 
cannot be explained by this phenomenon. The modified 
Kawaler braking law developed in response to problems 
in reproducing empirical trends, and a com parable mod- 
ification of the IReiners fc Mohantvl (|2012l ) prescription 
could be viable and interesting to explore. 

3.3. Amplitudes and Age Uncertainty 

As demonstrated in the previous section, a regime ex- 
ists where the range of stellar rotation rates has con- 
verged enough to make gyrochronology a good clock. 
For stars of the same mass, the age constraints from gy- 
rochronology are tighter for more slowly rotating stars. 
Because magnetic activity, and therefore the presence of 
starspots, declines with time, so too does the amplitude 
of photometric variability. This raises the question of 
whether there exists a good mass window where stel- 
lar amplitudes of variability are high enough that cur- 
rent ground based surveys can obtain rotation periods 
for the older, less magnetically active stars that are the 
most promising candidates for gyrochronology. In the 
following discussion, we demonstrate that current ground 
based studies can measure rotation-based ages for stars 
up to 1-2 Gyr. 

From the ground, current rotational studies have 
a minimum detectable amplitude of 3 m illimag (e.g. 
lHartman et al.ll2009t iMeibom et~aT1l2011al) . which may 
be pushed to as good as 1 millimag in the future. 
The space-based Kepler mission is monitoring roughly 
150,000 stars with a time cadence of 30 min. Its photo- 
metric precision is estimated to be 1 millimag for a 15.5 
mag star and as good as 10~ 5 mag for an 8 th mag stafl 
Indeed, some studies of ste llar amplitude from th e first 
quarter of Kepler data (e.g. iBasri et aLll2010l . 120111) have 
already begun examining how the amplitude of variabil- 
ity depends on T e g, periodicity, and plan to investigate 
the age distribution of these stars once a longer time 
baseline is available. 

In light of these upcoming developments, we use the 
current observations from open clusters to estimate what 
trends should be visible in the Kepler sample. To do so, 
we requi re a conversion between rotation period and am- 
plitude. lHartman et al.l (2009) used the clean sample of 
rotation periods in M37 to explore the relationship be- 
tween observed amplitude of variation and Rossby num- 
ber. Dividing their sample at (B-V)o = 1.36, the blue 
stars display a strong anti-correlation between Rossby 
number (i?o) and amplitude which flattens out for Rq < 
0.3. Redward of that division, the anti-correlation is less 
pronounced and the constraints on the convective over- 
turn timescale are weaker. Focusing on the bluer stars, 
lHartman et al.l call into question the period measure- 
ments for stars with Rq > 0.6, raising the possibility 

2 http://keplergo.arc.nasa.gov/CalibrationSN.shtml 
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Fig. 15. — Top: Maximum observable rotation period for a sur- 
vey capable of detecting amplitude variations of 10 (black; solid), 
3 (black; short dash), 1 (gold; dot), and 0.1 (red; long-dash) mil- 
limag. Bottom: The associate median age of a star with maximum 
period given in the top panel. Each line represents the maximum 
age a survey of a given precision could expect to obtain via gy- 
rochronology. Arrows indicate the direction of the trends. The 550 
Myr age bound is shown as the horizontal black dashed line. The 
HATNet survey (black bar) is capable of detecting the rotation pe- 
riods of 20 days at amplitudes of 0.01 mag, which corresponds to 
ages of 2 Gyr or younger. 



that they reflect spot evolution ti mescales rather than 
rotation periods. Based on this, IHartman et al.l find 
an empirical fit to the r-band amplitude of M37 cluster 
members with (B-V)o < 1-36 and Rq < 0.6 of the form 



A r = 



0.078 ±0.008 



1 



-Ro 



0.31±0.02 



3.5±0.5 ' 



(9) 



We invert Eq. [9] to determine the maximum observ- 
able rotation period for a survey of a given photomet- 
ric precision. For some amplitude A r , the resulting R 
may be converted into a rotation period if the tcz is 
known. We perform this calculation for a grid of dif- 
ferent stellar masses where we utilize the tcz predicted 
by the YREC models. Because tqz changes only very 
slightly over a star's main sequence lifetime, we adopt a 
mass- Tcz relation at 1 Gyr. We restrict the grid to the 
range where the Hartman et al. relationship is applica- 
ble. To determine this, (B-V) is converted to mass using 
a 1 Gyr, [Fe/H] M37 = +0.045 YREC isochrone. Un- 
der this transformation, the color limit (B-V)o < 1-36 
corresponds to masses M> 0.65 M©. This results in a 
mass-period relationship at fixed amplitude. These iso- 
amplitude curves represent the longest period for which 
a star of that mass could be detected with that photo- 
metric precision. A more slowly rotating star tends to 
have smaller variability in accordance with the activity- 
rotation-age relationship and therefore will drop below a 
given detection threshold. 

These lines of maximum period are shown in the top 
panel of Figure [__]_] Because older stars rotate more 
slowly, if a survey is unable to detect periods longer than 
some threshold, this places an upper limit on the age 
of stars that can be determined with gyrochronology for 



that survey. We translate the maximum rotation periods 
into age using the median of the rotation distribution at 
a grid of ages given in Table [2j The mass range of these 
calculations is set by where Eq. [9] is applicable and by 
the scarcity of high mass stars in M37. 

These maximum ages, shown in the bottom panel of 
Figure 1151 indicate that the current photometric pre- 
cision of ground-based surveys (3 millimag) limits gy- 
rochronology to measuring stars younger than ^1 — 2 
Gyr. An improvement of a factor of 3 in precision (~ 1 
millimag) would open up the possibility of measuring 
rotation-based ages for stars as old as 2 — 3 Gyr. A 
further order of magnitude improvement in precision is 
required to use gyrochronology to measure the ages of 
stars over the full age range of the Galactic disk (up to 
10 Gyr). Kepler is capable of achieving the requisite 
precision to do so for stars brighter than niKcp ~H mag 
ijBasriet al.ll20Tot) . 

The Ro- amplitude relationship suffers from a very lim- 
ited dynamic range in Ro from the M37 data. New ob- 
servations of field stars will extend the relationship to 
larger Rq (longer rotation periods) and provide a bet- 
ter constraint on the unsaturated slope. For now, we 
motivate a relationship with a shallower slope and con- 
sider its ramifi cations. This analysis is based on the 
IHartman et al.l best fit curve that maps Ro to ampli- 
tude. Some of the scatter in the M37 data is stochastic. 
The amplitude of variability of a star depends not only 
on the spot covering fraction, but also on the degree to 
which the spots are isotropically distributed over the sur- 
face of the star. At one extreme, if a star had a perfectly 
isotropic distribution of spots, no modulation in the light 
curve would be observed. Alternatively, large amplitude 
modulations would arise in the case of a concentration of 
spots along one line of longitude. By this logic, the up- 
per envelo pe of the amplitude- Rp diagram for M37 (see 
Figure 17, Har tman et al .1120091 ) represents the maximum 
possible amplitude of variability for a star of a given Rq. 
The cloud of points under this envelope could be pro- 
duced by different spot configurations that decrease the 
amplitude of variation. If instead we had used the up- 
per envelope, rather than the best fit through the cloud 
of points, to define the the transformation from Ro to 
amplitude, visual inspection indicates that the slope of 
the anti-correlation would be shallower. If this is the 
case, then using the best-fit -Ro-amphtude relationship 
would represent an underestimate of the amplitude of 
variability. That is, there may be older, more slowly ro- 
tating stars which exhibit higher amplitude variability 
than one would expect from Eq. [i__| because they pos- 
sess a more favorable configuration of spots. This would 
enable ground-based surveys with sufficiently long time 
baselines to probe older ages than predicted in Figure 
1151 Conversely, stellar variability may be suppressed to 
arbitrarily low levels by an nearly isotropic distribution 
of spots and therefore one should not expect all low am- 
plitude stars to be old. 

We can perform an empirical check on these maximum 
perio d-age-sensitivity predic tions using current field star 
data. Hartma n et al.l (J2011D searched for stellar variabil- 
ity in light curves collected by the HATNet survey for 
transiting extrasolar planet. This yielded rotation peri- 
ods for 1568 field F, G, K, and M dwarfs. HATNet offers 
a photometric precision of a few millimag for Re ~ 8 
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mag and 0.01 mag fo r Re ~ 12 mag. Figure 18 in 
lHartman et al.l (|2011f) compares the mass-rotation distri- 
bution for field stars detected with HATNet with those 
from open clusters, including M37. At 0.7 M Q , the field 
stars have rotation periods a factor of ~ 2 longer than 
the correspond ing stars in M37. This also agrees with 
lHartman et all s Figure 17 that shows that the number 
of detected stars rises slowly as a function of rotation pe- 
riod, peaks around 20 days, and then drops off steeply to 
larger periods for all color bins. If HATNet can detect ro- 
tation periods of ~ 20 days or shorter, this corresponds to 
ages K-star ages 2 Gyr or younger from interpolation on 
Table[5J We overplot the HATNet detection capabilities 
in Figure [15] and discover that they more closely match 
the curve at a mag = 0.003 mag than the 0.01 mag curve 
that corresponds to the surveys detection threshold. This 
is an encouraging sign that our predicted detection limits 
are actually too pessimistic and the activity-rotation-age 
relationship does not fall off as steeply as Equation |H] 
predicts. 

The results from HATNet mean that current ground- 
based surveys can identify populations of stars younger 
than 2 Gyr based on their amplitude of varability alone. 
However, such a sample could suffer contamination from 
merger products or synchronized binaries that are more 
active than isolated field stars of the same age. Addi- 
tionally, the detection efficiency among slowly rotating 
is unknown. Stars inclined pole-on to our line of sight 
or those at the minimum in their stellar activity cycle 
contribute to a non-detection fraction of stars at a given 
period. The impact of these stellar population effects 
will be discussed in 52) 

4. DISCUSSION AND CONCLUSIONS 

Our goal is to quantify the strengths and weaknesses 
of rotation as a stellar age indicator. For which stars is 
rotation a good chronometer, and how reliable are those 
ages likely to be? How does it compare to other age di- 
agnostics? Rotation decays as stars age, and old stars 
are therefore harder to detect than young stars. How 
does the age limit for photometric surveys depend on the 
precision of the photometry? In addition to these ques- 
tions, there are broader issues that deserve comment. 
Our data is concentrated in young stars, while the most 
promising applications are for old stars. What potential 
issues could complicate the interpretation of rotation in 
older stars? How can observational selection effects be 
accounted for? And what backgrounds need to be un- 
derstood for stellar population studies? In our summary 
we begin with the quantifiable answers to the first class 
of issues, and then turn our attention to the challenges 
and opportunities for progress in the future. 

4.1. Current Results 

Our main result is that stellar rotation is a promising 
tool for estimating the ages of Sun-like stars. The upper 
mass limit for rotation-age relations is set by the F star 
transition from deep to shallow surface convection zones. 
It is not clear that stars more massive than M ~ 1.2 
Mq have a rotation-age relationship, and it is difficult 
to detect rotation periods signals for stars more massive 
than the Sun because the fractional star spot coverage, 
and the amplitude of variability, drops steeply for hotter 
stars. 



We demonstrated the method of deriving rotation- 
based ages from spin down models works for stars below 
the F star transition and above the fully convective limit. 
Although both theoretical angular momentum loss mod- 
els and the semi-empirical published gyrochronology re- 
lationships are calibrated to the same open cluster data, 
systematic errors arise in the extrapolation of these re- 
lationships forward in time. Because of the limited data 
at late times, constraining the rotational behavior of old 
stars is difficult. We find that both a mod ified Kawaler 
law and the Mamajek & Hillenbrand (2008) gyrochronol- 
ogy relation agree with constraints from old field binaries, 
although their predicted ages at fixed rotation period dif- 
fer by as much as 30%. We view this as a measure of the 
current systematic error level in rotation-age relations. 

These loose constraints from old binary stars opens 
the door for other angular momentum loss laws offer- 
ing a different spin down trajectory than the modified 
Kawaler law. A reassessment of the 24 year old Kawaler 
braking law i n the modern context is wel come. However, 
the proposed iReiners k, Mohantvl (120121) i s incon sistent 
with data. Taking the IReiners fc Mohantvl (|2012l ) law at 
face value indicates that there are a strong mass- and 
radius-dependent trends in the spin down behavior at 
fixed w cr it. The torque is predicted to be so inefficient 
that it conflicts with long period rotation measurements 
from the old binary star systems (61 C yg and a Per) , an d 
field K a nd M dwarfs from CoRoT (JAffer et al.ll2012h . 
MEarth ( Irwin et all 120111) . and Mount Wilson survey 
(e.g. iDonahue et al.lll996fl. Future p hysically motivated 
revisions to the Reiners fc Mohantvl law may solve these 
problems. TheJ3arnesj (2007]) relation is also inconsistent 
with the binary star data. More rotation and age data 
is necessary for old stars to provide better guidance to 
theory. 

The age errors for gyrochronology relationships are 
larger than optimistic early claims. We investigated the 
impact of the range of stellar birth rates and surface lat- 
itudinal differential rotation on the precision of age mea- 
surements. These two phenomena will induce noise in 
age estimates of coeval stars even given perfect knowl- 
edge of the mean age-rotation relationship. Below 0.6 
- 0.7 M , the minimum errors in rotational ages rise 
steeply; a range of surface rotation rates persists even 
for old stars. The exact boundary depends on how the 
period data is transformed from the observational (color, 
period) domain into the theoretical (temperature, mass, 
age) domain. Fortunately, there is good empirical guid- 
ance for both the time evolution of rotation and the 
proper choice of color-temperature relationships for low 
mass stars. There is thus a window in the 0.6 - 1.2 M Q 
domain where stellar rotation is a promising chronome- 
ter. 

Our technique for assigning errors for ages derived from 
rotation is conservative, and it may be possible to ex- 
tract useful statistical age information from rotation in 
M dwarfs. The first step would be to empirically mea- 
sure the shape of the rotation distribution and to em- 
pirically calibrate the time dependence of the spin down 
in systems of known age, such as open clusters. One 
could then construct a probability distribution of ages 
for a fixed rotation period. In any given field, however, 
the true probability distribution would be a function of 
the age distribution in the background. A good exam- 



20 



Epstein & Pinsonneault 



pie would be a field with a star cluster superposed on a 
sparse background. A blind method would assign a range 
of ages to an intrinsically peaked age distribution. The 
age distribution in a given field would therefore need to 
be obtained independently (for example, with rotation- 
based ages of K stars) and a star by star age estimate 
could then be constructed for a sample of M dwarfs for 
which age estimates are otherwise difficult. This method 
would also be particularly useful for determining the ages 
of special targets, like planet hosts. 

Stars like the Sun do not have a unique surface rota- 
tion rate, and this injects another source of uncertainty 
in rotational age. The latter effect becomes dominant 
in old field stars, imposing a minimum uncertainty of 
20% in age estimates for a 10% range in surface rota- 
tion period. Different asymptotic angular momentum 
loss prescriptions have a major impact on the limiting 
precision, with the recent iReiners fe Mohantvl prescrip- 
tion implying much larger age uncertainties t han pre- 
vious spin down models based on a modified iKawalerl 
(1988) loss law. A proper empirical calibration of the 
late-time dependence should address the functional form 
of the rotation-age relationship, and it will also serve as a 
check on the width of the observed distribution. A proper 
calculation of the impact of surface differential rotation 
on period measurements would involve modeling of the 
number and latitude distribution of spots as a function 
of starspot cycle, and our method may overestimate the 
practical impact of surface differential rotation. A larger 
problem may be the poor visibility of solar analogs during 
minima in their activity cycles (see below.) Because the 
timescale for stellar activity cycles is of order decades, 
monitoring over many seasons (e.g. the Mount Wilson 
survey) would be necessary to measure a stars' full range 
in surface rotation periods. Nonetheless, surface differ- 
ential is not a fundamental limitation in the same way 
as the initial conditions are. With extremely long term 
monitoring, surface differential rotation could in princi- 
ple be calibrated out. In such a case, one might define a 
equatorial rotation period-mass-age relationship. How- 
ever, the time investment to do this for a large number 
of calibrating stars is significant. 

We also examined prospects for using rotational vari- 
ability to detect young populations from ground-based 
surveys. Current large ground-based surveys can reach 
a precision of 0.01 mag, and targeted open cluster cam- 
paigns have demonstrated 0.003 mag precision levels. We 
find that field K star populations younger than 1 and 
2 Gyr (respectively) can be detected at these levels of 
precision if we adopt the amplitu de-Rossby scaling re- 
lationship of Hart man et al.l (|2009f ). The field star data 
of Hartm an et al.l (| 2011| ) detected longer periods than 
would be predicted based on this estimate, suggesting 
that ground-based surveys may be capable of finding 
older stars than we predicted (2 Gyr at 0.7 M Q and 0.01 
mag precision, as opposed to our 1 Gyr estimate) . Much 
higher levels of precision, currently only possible with 
satellite data, are required to detect rotational modu- 
lation in old stars. With new ground-based surveys on 
the horizon, such as LSST and the Las Cumbres Global 
Telescope Network, the monitoring large fractions of the 
sky will reveal variability for large samples of field stars. 
Because the amplitude of variability, like rotation rate, 
declines with age, such surveys will be able to identify all 



young K stars by variability alone. A subset of active G 
stars will also be discovered, consisting of a very young 
(< 1 Gyr) sample and tidally synchronized binaries. 

4.2. Cautions and Future Prospects 

Our method implicitly assumes that we can safely 
project the observational pattern seen in young sys- 
tems into the future. There are interesting theoretical 
classes of models where this is not the case. For exam- 
ple, stars could vary in their core-envelope coupling be- 
cause of differences in internal magn etic field morphology 
(jCharbonneau fc MacGregorlll993f ). This would lead to 
a persistent range in rotation rates at late ages gener- 
ated by star-to-star differences in core rotation and the 
effective moment of inertia subject to spin down. It may 
not be possible to maintain a dynamo if the rotation 
rate beco mes too long relative to t he convective overturn 
timescale (jDurnev fc Latourlll978[ ): both spin down and 
starspots could simply cease in stars above some critical 
age for a given mass. Finally, gravitational settling could 
induce a composition gradient sufficient to decouple the 
surface convection zone from the i nterior. There is 
some evidence from old open clusters ([Sestito fc Randichl 
2005) that the pattern of lithium depletion stalls for sys- 
tems older than a Gyr; if rotational mixing is ti ed to 
angular momentum loss (|Pinsonneault et al.l 119981 ). this 
could be evidence that there is a similar phenomenon in 
the torque. These questions can be addressed with Ke- 
pler data in old open clusters and targets with known 
seismic ages and measured rotation rates. 

Interpreting rotation in field stars will also require care- 
ful consideration of observational selection effects and 
stellar population backgrounds. It will be difficult to de- 
tect stars observed nearly pole-on (Krishnamurt hiet al.l 
119981 : lHartman et al.l [2010 : Uackson fc Jeffries! 120101 ) . It 
will also be difficult to detect solar-like stars either at 
the minimum of their sunspot cycle or if they are in ex- 
tended Maunder minimum-like epochs. What is even 
more challenging will be modeling the stochastic effects 
of the distribution of star spots on the surface and obser- 
vational selection effects from the observing cadence. All 
of these effects conspire to produce a substantial popula- 
tion of non-detections, and careful work will need to be 
done to quantify their impact. 

Finally, our age estimates, as well as the body of gy- 
rochronology work, assumes that we have a sample con- 
sisting of main sequence single stars. There will be sig- 
nificant background populations of tidally synchronized 
binaries and blue stragglers. A few perc ent of solar-type 
stars are likely to be merger products (|Andronov et al.l 
|2006l ), and their rotation rates will not be simple func- 
tions of their age. Subgiants may be difficult to distin- 
guish from dwarfs, especially for F and G stars, and their 
rotation-age relationships will be very different from one 
another. All of these backgrounds will produce spuri- 
ous signals when mapped onto single star relationships. 
Finally, the metallicity dependence of spin down rela- 
tionships is poorly understood because the calibrating 
open cluster sample has a very limited dynamic range in 
composition. 

4.3. Gyrochronology in Context 

Using rotation as a clock is useful if it is competitive 
with other techniques for determining ages for unevolved 
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Fig. 16. — Comparison of expected age precision for a 5 Gyr 
old star from different techniques of estimating the age of main 
sequence stars below the turnoff. The contenders include rotation 
based ages (solid), asteroseismology (dashed), and trigonometric 
parallax combined with [Fe/H] measurements good to 0.1 dex (dot- 
ted) and 0.05 dex (dot-dash). 

dwarfs. Ages may be determined by comparing a star's 
location on the CMD to stellar interior model predictions 
if the distance and metal content of a star is known. Gaia 
will provide exquisite geometric parallax measurements 
accurate to a few /xas for ~ 10 9 stars. However, with- 
out helium or mctallicity information, stellar isochrones 
suffer a degeneracy with age. Figure H6l shows the possi- 
ble age precision achievable with [Fe/H] good to 0.1 and 
0.05 dex, optimistically assuming perfect distance, tem- 
perature, and helium measurements. To generate these 



curves, we use Dartmouth isochrones to determine what 
difference in age produces a change in luminosity equiva- 
lent to a shift of ±0.1 dex (or 0.05 dex) in [Fe/H] at fixed 
temperature. Gaia's photometric accuracy will make the 
precision of metallicity information the limiting source of 
error. 

With the advent of Kepler, asteroseismology is expe- 
riencing a renaissance. Asteroseismology obtains age 
estimates for dwarfs by measuring the helium content. 
These arc good to 10% of a star's main sequence lifetime 
(Kicl dsen et al.l [2009f ). making the best age constraints 
possible for more massive stars. Figure [16] shows that 
asteroseismology potentially offers better age constraints 
than isochrone-based techniques in the low mass regime. 
However, stellar pulsa tions may not even be detectable 
in low mass stars (e.g. iChaplin et ai1l2011| ). 

Because both geometric parallax and asteroseismology 
measure ages based on changes on the nuclear timescale, 
both suffer higher uncertainties for low mass stars. In 
comparison, Figure [TBI demonstrates that rotation consis- 
tently offers ages good to ~ 1 Gyr for stars M> 0.55 M . 
Rotation-mass-age relationships are not competitive with 
other techniques for stars more massive than the Sun 
and the spread in rotation is large for low mass stars. 
Nonetheless, in an intermediate mass regime, rotation- 
based ages are not only competitive, but clearly the best 
possible age indicators available for field dwarfs. 

We are grateful to Donald Terndrup for his thoughtful 
comments and illuminating discussions in the prepara- 
tion of this manuscript. C.E. and M.P. were supported 
by NSF grant AST-0807308 to The Ohio State Univer- 
sity. 
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